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\ ABSTRACT 

We report the "state of the art" of the problem of B + L violation in high-energy 
1 electroweak scatterings. Results of various analyses point toward (though do not prove 

rigorously yet) the "half-suppression", i.e., that the B-\-L violating cross section remains 
suppressed at least by the negative exponent of the single instanton action, at all energies. 
^/■^ , Most interesting techniques developed in this field are reviewed. Particular attention is 

■ paid to unitarity constraints on the anomalous cross section, and to some conceptual 

problem involving the use of the optical theorem in the presence of instantons. 



1. Introduction 

Oh, We first present an overview of the progress made in the last few years in the 

problem of B + L violation in high energy electroweak processes. The earlier part 
of the development will then be reviewed, which serves as a technical introduction 
to the subsequent sections. 
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H , 1.1. Overview 



In the standard electroweak theory, the baryon (and lepton-) number is not 
strictly conservectl as a consequence of the chiral anomalyB 



d ^ = T^Tri^i^ (1) 



16tt 2 
where 

J" = ^Ll^L (2) 

is the chiral current for a lefthanded doublet ipL in the theory. Due to dynamical, 
non perturbative effects (such as instantons or sphalerons) this leads to physical 
processes with the selection rule, 

A(B + L) = (3 + 3)AA/bs; A(B-L)=0, (3) 



i 
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where Afcs is the Chcrn Simons number (see Eq.(178)). For instance, an instanton 
leads to elementary processes such as 

Ql + Ql^ 7q L + 3l L + n w W + n z Z + n h H. (4) 



The instanton describes a transition between the neighboring classical vacua 
It is well known that in the standard electroweak model the barrier height corre- 
sponds to the energy of the sphaleron configuration!, 

E sp = -M w ~ lOTeV, (5) 
a 

(where a = a\y ~ 1/32.) Therefore one expects that at energies much lower than 
E sp the cross sections for such A(B + L) ^ processes are typically suppressed by 
the t'Hooft's tunnelling factor, 

e -4f „ io-i7o (6) 

utterly too small for them to be observable. 

Nevertheless, in at least two situations, at high temperatures0 and at very high 
fermion densitiesEl, baryon number violation is believed to proceed unsuppressed. 
The interest and importance for such a possibility is mainly related to the problem 
of cosmological baryon number generation in the standard SU(3) x SU(2) x L/Xl ) 
theory of fundamental interactions. For a recent review, see e.g. Shaposhnikovca. 

Whether or not the A(B + L) ^ events will become observable at high energy 
scattering processes (at the energy available at SSC or LHC), has been the issue of 
an active debate for a few years now. From experimental physics point of view, the 
prospect of observing roughly isotropic .production of a large number n = 0(— ) of 
W or Z bosons is quite an exciting oneQ. From the theoretical point of view, the 
problem involves some of the most subtle aspects of non-Abelian gauge theories, 
such as chiral anomaly, degenerate classical vacua and quantum tunnelling among 
them, large order breakdown of perturbation theory, compatibility of semiclassical 
expansions with unitarity, and so on. 

It is the purpose of the present paper to review the latest developments in this 
field of research and assess our general understanding. Earlier works have been fully 
reviewed by Mattisil. See also the proceedings of the Santa Fe workshop!. For a 
more recent review, see Tinyakovliil. 

Historically, after first suggestions0 : 0, t he s em i-cl assi cal estimate of the total 
cross section with A(B + L) =/= was shown^LLito grow exponentially with energy 
(at low energies), stimulating further workscj til. It was noted immediately^ that 
the growth could not continue indefinitely as the computed cross section violated 
the unitarity limit above the sp hal ero n mass energy E sp . Quantum corrections 
around the instanton were studiedE3~EJ and it was shown that all tree type correc- 
tions involving the final states (so called soft-soft corrections) contributed to the 
nontrivial, exponential energy dependence. In particular, it was arguedtJ that these 
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corrections could be expressed in the form, 

<7A(B+L)^ = e^ F ^ (7) 



where the function F(x) (called sometimes "Holy Grail" function) has an expansion 
in x = E/E sp «1. 

Two distinct approaches, the R term methocO'EZI and the valley methocEj o, 
were developed to com pute the anomalous cross section, and used in recent per- 
turbative calculationsE3~El Subsequently the equivalence of the two methods as 
regards the final state corrections, has been argued to hold to all orders of pertur- 
bations around the instantorH. The developments which followed, to be discussed 
in Sections 2 and 4, make crucial use of these methods. 

A non perturbative model, utilizing the valley field and the optical theorem, 
was presented by Khoze and Ringwaldcj, with the aim of resumming all final tree 
corrections. 

Other kinds of corrections involving the initial, hard particles (the so called 
hard-hard and hard-soft corrections) were then studiedEj E3, following the impetus 
provided by the works of Mueller. The results suggest the exponentiation of loop 
contributions to all orders, however, in the direction of suppressing the A(£> + L) 
cross section. 

Such an exponentiation supports strongly the idea that the quantum corrections 
involving the initial hard particles can also be included in a modified semi-classical 
approximation. Led by this thought several different approaches have been pro- 
posed. Works reviewed in Section 2 may indeed be called "Search for the new 
semi-classical field" , or perhaps, more poetically, IlSearch for the Holy Grail." 

MuellerE^I and independently Mc Lerran et al.ESl, propose a classical equation of 
motion for the fields, with the source term to take account of the initial energetic 
particles. They show that to lowest orders their solution automatically reproduces 
the quantum loop corrections found earlier by direct calculations. However, as 
pointed out by these authors themselves, the necessity of using Minkowskian time 
or complex time (and in general complex fields) and the consequent complexity of 
the field equation involved, seem to make the task rather a formidablaxineL 

An alternative approach was proposed by Rubakov and TinyakovlI3~E3. They 
propose to study, instead of the original 2 — > all type cross sections, the processes 
involving initial many-body states, n — ► all, n = -3, in the limit g 2 — > 0, with fixed v 
and fixed total initial energy E. The idea is that once a semi-classical approximation 
for such many-body processes has been established, one can then take the limit 
v — > 0, hopefully recovering the answer for the 2 — > all process. The coherent state 
formalism first introduced in the A(B + L) =^ problem by Khlebnikov and others^ 
plays a powerful role in this approach. Latest workE.3 seems to show that the above 
limit is indeed a smooth one. 

It is however not clear at the moment whether these approaches can lead to a 
new truly semi-classical approximation and what the eventual answer might be. 

Perhaps the clearest physical picture of what happens in high-energy A(B + 
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L) 7^ Dip cesses, has been given by the study of quantum-mechanical analogue 
problemsEj E3. There are two competing factors which determine the energy de- 
pendence of the cross sectionnJ. One, the semi-classical tunneling factor, rapidly 
grows with energy, and the suppression is altogether lifted when the sphaleron en- 
ergy is reached. Another factor is the probability amplitude - mismatch - between 
the most favorable state (for the purpose of tunnelling) and the initial two-particle 
state. And this factor (an analogue of Landau's semiclassical matrix elements) 
gets strongly damped as the energy increases. As a consequence the full amplitude 
never gets large, leading to the "half-suppression" result (see below). Diakonov 
and PetrovH recently applied to the problem a generalization^! of such a WKB 
approximation to field theory, finding some indication for the suppression of the 
cross section for the isotropic production of many VF's, even at very high energies: 
E 3> E sp . 

Quite parallel to the developments mentioned above, several argumen ts, es- 



sentially all taking unitarity costraints into account, were presentedEEl Ei3, which 



suggested that the A(B + L) ^ cross section remains suppressed at least by 

e-S (8) 

(so-called "half suppression"). In spite of their apparent differences, these argu- 
ments are actually closely interrelated, and depend only on the S-wave unitarity 
and the dominance of multiparticle production events. 

The who le question was analysed from a somewhat different angle by the present 
authorsEi3'E23. The approach using the optical theorem and the valley field to com- 



pute the total anomalous cross sectioned, takes unitarity into account automatically, 
and furthermore seems to enable one to compute the Holy Grail function up to the 
sphaleron energy. However, a new kind of problem arises (which may be termed the 
"unitarity puzzle"). Namely, how can one extract the part of the imaginary part of 
the elastic amplitude, that corresponds to the A(B + L) ^ intermediate states? 
In other words, how is unitarity satisfied in the presence of topologically nontrivial 
effects such as instantons ? A partial answer was given by the equivalence proof 
by Arnold and MattisEl; however the crucial issue concerns the initial particles, see 
Section 5. Or, when does an " instanton-anti-instanton" type configuration cease to 
be topologically nontrivial? To clarify these issues requires a detailed study of the 
behavior of chiral fermions in a background of instanton anti-instanton type (hence 
the title of this review!). The result of this investigation leads us once more to the 
above mentioned "half-suppression" of the A(B + L) ^ cross section. 



The rest of the paper is organized as follows. In subsections 1.2, 1.3 and 1.4 
the earlier developments, including the original instanton calculation, the i?-term 
method and the valley method, are reviewed, which also serve for fixing the con- 
vention and notation for later sections. The reader already involved in the research 
in this field may well skip this section; for others this section should provide an 
appropriate technical introduction. 
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The Section 2 contains a somewhat detailed review of more recent, new semi- 
classical approaches to the Holy Grail function. The discussion is divided into 
two parts, the first (subsection 2.1) dealing with attempts to take into account 
the effects of the initial high energy particles into a semi-classical equation, the 
second (subsection 2.2) being mainly concerned with the multi-particle approach of 
Rubakov and Tinyakov. 

In Section 3 we review those studies based on quantum mechanical analogue 
problems, which appear to provide an intuitive understanding of the whole problem. 

In Section 4 various analyses, leading to unitarity bounds on the baryon number 
violation in high energy electroweak scatterings, are reviewed. We first discuss 
a simple, multi-instanton unitarization picture and the general unitarity bound 
following from the S-wave dominance (subsection 4.1). The results are corroborated 
by an explicit resummation of multi-instanton contributions in subsection 4.2. The 
physical difference between the high temperature or high density transitions and the 
high energy B + L violation is briefly mentioned in subsection 4.3. 

The calculation of the anomalous cross section via optical theorem, done with 
the valley method, is critically analysed in Section 5. After the discussion of the 
unitarity puzzle (subsection 5.1), the fermion Green function in the background of 
the valley field is studied and the anomalous part of the forward elastic amplitude 
identified (subsection 5.2). The melting of the instanton anti-instanton pair, and the 
ensuing transition to purely perturbative amplitude, is discussed in subsections 5.3 
(through the study of the Chern-Simons number) and in 5.4 (in which the fermion 
level crossing in the valley is analysed). The implication of these results to high 
energy electroweak processes is summarized in subsection 5.5. 

We conclude (Section 6) by discussing a unifying and consistent picture which 
seems to emerge through different types of analyses reviewed here. 

1.2. The original instanton calculation 

In this section we recall briefly the calculation of RingwaldE] and EspinosaB for 
the cross section of the process (0). See MattisEl for more details. In order to keep 
formulas as simple as possible we shall take only fermions and W bosons as external 
particles. Let us consider the Euclidean n + Np point Green function (Np being 
the number of the lefthanded fermion doublets), 

G(xi, . ..x n) yi, . ..y NF ) = 

J dAd^ d^d0A^( Xl ) . . . A^x n )^( yi ) . . .i>(y NF ) e~ s ^ A ^^ . (9) 

This Green function gets contributions only from gauge fields with unit topological 
number Q = 1, where 

Q = f d 4 x-^Tr F^Fp, = M CS (~oo) - M CS (oo), (10) 
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(Afcs is the Chern Simons number, see Eq.(178)). In the semi-classical approxima- 
tion one needs appropriate classical solutions with a finite action. In the elcctroweak 
theory one can show that no classical solutions exist in the presence of the Higgs 
vacuum expectation value. Indeed if one scales 

4>{x) -> 4>(ax); A^x) -> aA^ax), (11) 

(which is a particular kind of variation), and considers the action as a function of a, 
one sees that no minimum exists for a ^ 0, for any functions <fi and A. On the other 
hand the existence of a classical solution would require the action to be a minimum 
under all variations around it. 

One is forced to look for a constrained instanton solutionEl. Although the solu- 
tion is not known in a closed form, one knows asymptotic forms when the distance 
from the center of the instanton Xi is much less or greater than the size p of the 
instanton. Very near the instanton center it is given by 



g v p " {x-xif{(x-xi)* + f?y 

i / /'.-'/ ! , ,..,-1 \ / {x-Xj) 2 \ ^ (x - Xj) 2 

* {x) ^ u[ v/v2){((x-x i r + ^) ' ^^ <<:1 ' (12) 



where = (i, <?), = (cr^ and U is the global SU(2) rotation. Far away 
from the center A]^ nst ^ (x) is proportional to a massive boson propagator; 0( ms *) (x) 

(iUSt) ' r-r I 



approaches a constant, cf> (mst >{x) ~ U \ ^^y^ 

The leading order approximation consists in keeping only up to the quadratic 
part of the fluctuations in the action and in substituting the fields by the classical 
ones in the pre-exponential factors. For fermions the "classical solutions" are the 
zero modes of the Dirac operator in the classical background, which behave as a 
free massless propagator asymptotically. The gauge field satisfies the massive free 
field equation far from the instanton center: its Fourier transform displays a pole 
at k 2 = My/- Through the LSZ procedure, one gets the residue of this pole equal 
to 

fc-(p,k) = — P %„K, (13) 

where fj is the usual t'Hooft symbolR Each fermion doublet contributes with a 
factor ~ e lkxi p, coming from the LSZ amputation of the zero mode, 

ipotx)^-. — for x — > oo. (14) 

(a; - XiY 

The action evaluated at the classical field, is equal toll 

S c = — + 7t 2 p 2 v 2 + 0{p A v 2 M 2 H \og{M H p)). (15) 
a 
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Putting all these pieces together, and integrating over the collective coordinates 
(the instanton position, size and the isospin orientation U) one finds 

A 2 ^ n ~ ex P (-S mst ) | kx | .... | k„ +2 | (16) 

v 9 

where Si nst = 2ir/a and n\ comes from the integration over the size, 

2n 



dpp zn e-* " v . (17) 

This gives the S matrix element at fixed n. The total cross section is found 
by squaring it and summing over the final states. To understand qualitatively its 
behaviorS, recall that the relativistic phase space goes as E 2n /n\; the squared 
amplitude as (n\) 2 E 2n , where E ~ E/n is the average momentum. For large n 
these together_giiLe rise to an "exponential" series X^n^ 4 ^ 3 ) 3 ™ / (3ft)!. A more careful 
estimatioi£3'ESE3 leads to 

1 V > 1 y 3S —2Si ns t 

CT 2^all « ^p3i g j2 T 2^) e 

„ i e *(-i+K*)*), (18) 
s 

where i?o = VGitMw fee- 

The cross-section grows exponentially with energy. If such a growth should 
continue up to and above the sphaleron energy, 

E sp ee -M w , (19) 
a 

the exponential suppression factor would be compensated altogether. Actually, 
the approximations leading to Eq.(|l^) are valid for energies much less than the 
sphaleron energy as will be seen below. 

A simple calculation shows that the cross section is dominated by the production 
of W bosons, with the average number 

n^-x^, (20) 
a 

where the definition 

x=^~ (21) 

will be used throughout this review, x here must be small for these approximate 
estimates to be valid, but should not be too small either, so that the average mul- 
tiplicity ( pp| ) is sufficiently larger than unity. Eq.(^0|) then implies that the mean 
energy of the final particles is E/n ~ Mn/i -1 ' 3 , showing that the final state parti- 
cles carry soft momenta, while those in the initial ones hard (large) momenta. 

The inclusion of the Higgs particles does not change these features in any essen- 
tial manner. 
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Eq.(18), if extrapolated above the sphaleron energy, violates the unitarity limit. 
It must therefore be substantially corrected before the initial energy reaches E sp . 
Does the A(B + L) ^ cross section nonetheless grow up to the point of saturating 
the unitarity limit, i.e., to a geometrical size, ~ 1/s? We are interested here in those 
corrections which exponentiate and modify the energy dependence in the exponent 
of the result, i.e. , those which contribute to 

F(x) = Umlog^ajj. (22) 

F is called the Holy Grail function. 

The corrections to the leading instanton result, are often classified into three 
groups (see Fig. 1): "soft-soft", "hard-hard" and "soft-hard" corrections. The first 
type involves only the final particles, the second one the initial particles and the 
third one both the initial and final particles. Also, contributions coming from other 
classical solutions with unit topological number such as multi-instanton configura- 
tions, could be important (see Section 4). 

All the tree soft-soft corrections have been shownS to exponentiate and give rise 
to 1/q terms in the exponent. a In order to compute them, two distinct methods have 
been used: the valley method and the "i?-term" method, which will be reviewed 
briefly in the next subsections. 

1.3. The R-term method 

In the i?-term method, developed by Khlebnikov, Tinyakov and RubakovEl, the 
summation over the final states is taken into account as a sort of correction to 
the action. This could be convenient for an application of the steepest descent 
methodES. We follow here an approachto the method developed by Mc LerranE3, 
which is simpler then the original oneEj'El The cross section is expressed as a double 
path integral, which within the perturbation theory is nothing but a compact way 
of implementing the standard Cutkovsky's rules. 

Consider a theory with a scalar field <f> with Minkowskian action 

S = J d*x{ l -(d^Y l -m 2 4> 2 - V (</>)} (23) 

A Green's function (in momentum space) containing two initial particles and n final 
particles can be expressed as a path integral 

G(pi,p 2 ,kx,...k n ) = j dcf><f>(p 1 )*<f>(p 2 )*ct>(k 1 )....cj>(k n )e iS W. (24) 

To get the total cross section one applies the LSZ procedure to the above, squares 
it (by doubling the fields) and sums over all possible final states. These operations 



a For x = O(l) the naive perturbation theory around the instanton would give an increasingly 
divergent series of type ~ [ — ) n which should necessarily be summed to all orders, i.e., nontrivial 
corrections to the Holy Grail function. 
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are compactly expressedE-j by the introduction of 

R(k) = lim (fc 2 -m 2 ) 2 (f)(k)<P*'(k). (25) 

k 2 — >m 2 

a tot c & n then be written as a double path integral: 

atot = i /# / d4>'R*( Pl )R*(p 2 )e iSW - iS ^ + f (26) 

where F is the flux factor, F ~ s. Note that, thanks to the statistical factor 1/nl for 
identical particles, the factors associated with the final particles nicely exponentiate 
upon summation over n, and give rise to the third term in the exponent, the so- 
called R - term. 

The i?-term method is particularly powerful, when applied to the computation 
of the corrections to the leading semi-classical result, E q. (|l8| ) . We shall here content 
ourselves however to observe that to lowest order Eq.(|26|) reproduces Eq.(|Ts|). 

In this approximation one splits the fields into the classical and quantum parts, 

<j>(x) =^ mst \x;a)+^(x), (27) 

where a indicates the ensemble of the collective coordinates, and substitutes the 
classical solution into the R - term, as well as in the ordinary actions and in the 
pre-exponents in (^6|) . In the case of the standard model, the leading contribution 
to the exponent comes from W bosons and is represented in Fig. 2. Using Eq.([l3|) 
and by carefully continuing from Minkowskian to Euclidean space, one gets 



otot oc^| < n|A|p,-p > | 2 = e ~ W7r2/92 f dt dx dp dp 1 C e w , 
n J 



(28) 



where (t, x) is the Minkowski continued difference of the instanton positions, Xi — x\ 
and C contains the functional determinant over the nonzero modes, the Jacobian 
associated to the introduction of the collective coordinates as well as factors from 
fermion zero modesLi 

The crucial exponent W is given by 



W 



-iEt ~ 7T 2 v 2 (p 2 + p' 2 ) + ^p 2 p' 2 [ —e 1 ^^ 2 + k 2 ), (29) 

g J uj k 



where u>k = "v/k 2 . 

Integration over x by the saddle point method sets x = 0, which leads to 

W = - t Et-n 2 v 2 (p 2 +p' 2 )+ m7T 2 P . P . (30) 

9 t 

The last integration over t and p, p' by the saddle point approximation yields 
the well known result, Eq. (|l8|) . 

Higher-order corrections can be systematically taken into account by using the 
expansion Eq.(p7|) in Eq.(p6|). The resulting series in the Holy Grail function turns 



10 Fermions, anomaly and unitarity 



out to be a power series in x 2 ^ 3 (see Eq.(|39|) below.) It means that in order to 
get a significant result for the energy of the order of E sp (x ~ 1), one must be 
able to resum all the R-term contributions corresponding to soft tree graphs, or 
equivalently to solve the complete classical equations, including also the R-term. 

It should be emphasized that the i?-term method (and the i?-term corrected 
classical equation) only makes sense in Minkowski space. For instance, by "instan- 
tons" one really means their analytic continuation to Minkowski space. Also, a new 
solution of i?-term corrected classical equation with an 0(3) symmetry ("distorted 
instantons") can be found (though only in the linearized approximation)^, but it 
turns out to be impossible to continue it to Euclidean space. As a result it is a 
nontrivial task to get a nonperturbative control of (even) the soft-soft part of the 
Holy Grail function (which we shall call F ss (x) to distinguish it from full F(x)). 

The i?-term method has been generalized by EspinosatlZl, so as to take account 
of the effects of fermion pair production. His calculation shows that the inclusion 
of the fermionic i?-term does not modify the result to exponential accuracy: it does 
not affect the Holy Grail function. 

1.4. The valley method 

An alternative app roach to the calculation of the total cross section is based on 
the optical theorernEEf'0 which relates it to the imaginary part of the forward elastic 
amplitude. However since we are here interested only in the inclusive A(B + L) ^ 
cross section, not really the total cross section, a highly nontrivial problem arises 
in extracting the "anomalous piece" from the full imaginary part of the forward 
elastic amplitude. This issue will be discussed extensively in Section 5. 

0"A(B+£)/o is given by 

VAiB+L)^ = -Anom Im LSZ Wick J dAdip dip ijjip'ip'ip e ~ s ^ A ^I ; (31) 

where with Anom, Im, LSZ , Wick we indicated symbolically the operations of ex- 
tracting the anomalous piece of the imaginary part of the on-shell amplitude (LSZ 
procedure), from the four point function Wick continued to Minkowskian space. One 
wishes to evaluate Eq.(|3l|) by a semi-classical approximation. Since one deals here 
with an elastic amplitude, the relevant gauge background must belong to the trivial 
sector, with zero topological number. At the same time, however, it must describe 
nonperturbative effects of producing A(B + L) 7^ intermediate states: it must 
be topologically nontrivial locally, as e.g., the instanton anti-instanton pair. As is 
well known from the example of a quantum mechanical double well, a simple sum 
of instanton and anti-instanton (at finite distances) is not an adequate background, 
since the expansion around it produces a large linear term in the "quantum" fluc- 
tuation, because is not a solution of field equations. The effort to minimize the 
latter term by a shift of the field introduces automatically an effective interaction 
between the instanton and anti-instanton. 

The valley (or streamline) method provides a way to take such interactions into 
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account systematically. The valley method0 _ is a generalization of the standard 
Faddeev-Popov procedure of treating zero modes in quantum mechanics and in 
quantum field theory, to the case of quasi zero modes. If the action has a valley-like 
shape in the functional space, i.e. if its value slowly changes along the bottom of 
the valley (streamline), one must first perform Gaussian integrations in directions 
orthogonal to the streamline, then integrate the result along the streamline. 
The valley trajectory (streamline) <f> a is a solution of the equation, 13 

w(s , a) __ = _|^ (32) 

(w is an arbitrary weight function) , a parametrizing the bottom of the valley. The 
generating functional Z in leading approximation is 



(33) 

This would be the same as the standard functional integration with a zero mode, 
were it not for the nontrivial integration over a. 

In pure Yang-Mills theory in four dimensions, the valley equation was solved0-@ 
by using an Ansatz, = (2/ g)d^ u {x u /x 2 )s{x 2 ). The problem is then reduced to 
that of a simple one-dimensional quantum mechanical double well, for which the 
solution to Eq.(^||) is known. The valley trajectory A^ allev ^ found this way (after 
a particular conformal transformation of the original ansatz) is 

i / _ r \ r V*^ ^a)f 



(x - Xj) v p 2 

(x - Xi) 2 ((x - Xi) 2 + p 2 ) 
(x-Xi+y) v {x-Xi) v 



(x-Xi+y) 2 (x-Xi) 2 



(34) 



We defined 



y" = -i?V(z-l); 



z 



(R 2 + 2p 2 + Vi? 4 + Vi? 2 )/2p 2 , (35) 



where x^ and are the centers of the instanton and anti-instanton, p is their 
(common) size. As is seen from Eq. (|34|) j\& alley ^ interpolates between two solutions 
of the classical Yang-Mills equation: the simple sum of instanton and anti-instanton 
at infinite separation (R — oo) and a gauge-equivalent of = (at R = 0). 

For simplicity of writing, the size of the instanton and that of the anti-instanton 
will be taken to be equal here; no generality is however lost since the saddle point 

'-'There is another jxwsible way to define the valley trajectory^ which is claimed to have some 
advantage over Eq. (p2|) . 
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equations for p a and p 1 set them equal anyway, in the problem one is interested 

here. 

The action of the valley is givenEj'EJ by c 



g ^{valley) \^ _ 



48?T 

~F 



6z 2 -U 17 f(z-5/z)(z+l/z) 



— — log z 



{z-l/z) 2 3 b V (z-l/z) 3 
which asymptotically (as R — > oo) behaves as 



(36) 



^ (37, 
a a ft 4 R° 

The first term is twice the instanton action, while the second one represents the 
attractive interaction between the instanton pair. 

The behavior of the valley action as a function of R/p is plotted in Fig. 3. 

In the case of the Weinberg-Salam theory, the solution of the valley equation 
is not k now n At energies much lower than the sphaleron mass, however, one can 
justifyS'BSthe use of Eq.@, Eq.@ in Eq.(|l]). Furthermore the fermion fields 
ip's or Tp's in Eq.(|3l|) may be replaced by the standard zero modes ip^'s (in the 
anti- instanton background) or t/Sq 's (in the instanton background), respectively, 
see subsection 5.2. In this manner the result Eq. (|l8|) , is reproduced by the valley 
method, confirming once more the exponential growth of cta(b+l)^o a ^ ^ ow ener gi es - 

The equivalence between the i?-term method and the valley approach has been 
derived perturbatively to all orders in x (or equivalently in J|) in a rather formal 
wayE3; their proof however neglects incoming particles (see Section 5) and also skips 
the problem of weight dependence of the result (see Section 2.1). 

Without knowing the solution of the complete valley equation, the utility of 
the valley method in the electroweak theory is limited, unfortunately. Nonetheless, 
Khoze and RingwaldE] attempted a nonperturbative calculation of the Holy Grail 
function by doing the following simplifications or assumptions (justified or not!). 
They (i) just add to the valley action Eq.(^6|) the Higgs contribution —2ir 2 p 2 v 2 ; 
(ii) substitute the fermion fields by the standard fermion zero modes ipQ^'s and 
"ipQ ^'s; and (Hi) evaluate the resulting integrations 

a = Im J dRdp cxp(ER - 2tt 2 P 2 v 2 - S vaUev (z)), (38) 

by the saddle point method (note that they commuted the Wick rotation with Im) . 
The saddle point equation relates the relevant values of p and R to the initial energy 
E. 

They find that the Holy Grail function increases monotonically and reaches 
precisely zero value (hence no exponential suppression of baryon number violation) 
at an energy of order of the sphaleron energy! (x = xkr = 8 a/3/5). See Fig. 4. 



c The fact that the action depends on the valley parameters only through 2 is a reflection of the 
conformal invariance of the classical Yang-Mills equation. 
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This result is repeatedly referred to in the literature as a very encouraging sign 
that the A(B + L) ^ cross section might reach the geometrical size at high 
energies: after all, if a reasonable dynamical model gives an interesting result, why 
could it not be true also in the real world? Unfortunately, there are strong reasons 
to suspect that the assumption (ii) used by Khoze and Ringwald is too naive and 
is hardly justified, precisely for values of R/ p < 1 where the valley action sharply 
drops to zero. We shall come back to the discussion of this model towards the end 
of Section 5. 

Note that, by construction, the valley field depends on the choice of the weight 
function w. other choices might possibly introduce different kind of difficulties 
such as the "bifurcation" (i.e. the loss of the saddle point of Eq.(0))EJ'E2|-Ell. A 
Minkowskian space formulation of the valley method was also givencJ. 

Finally l£t us mention the recent perturbative calculations of the Holy Grail 
function.E3~El Balitsky and SchafeiS use the valley method (they check the result 
by using the effective Lagrangian approach) to compute the third term of the Holy 
Grail function, with the result: 



where Eq = \f&E sp so that E/Eo — x/y/6. Silvestrov instead uses the R - term 
method to recover the coefficient of x 8 ^ 3 log(l/a;) , finding however only the piece 
(which agrees with Balitsky and SchaferEa) depending on the Higgs and W bo- 
son masses. Diakonov and PolyakovS get the other piece, but with a factor 1/2 
compared to Balitsky and Schafer (these authors however work in the pure gauge 
sector). 

2. In search of the Holy Grail. 

In this section several recent attempts to compute nonperturbatively the full 
Holy Grail function will be reviewed in some details. 

2.1. Initial corrections and attempts for a modified semi- classical ap- 
proximation 

Let us start with the discussion of corrections involving high-energy initial par- 
ticles, i.e. , corrections to Eq.(p"8|) due to interactions between hard particles or 
between hard and soft ones, (this definition can be applied indifferently both in the 
valley and i?-term approaches). The lowest hard- hard quantum correction in the 
simplified case of bosonic particles is shown in Fig. 5, and hard-soft (hard-hard) 
corrections in the fermionic case in Fig. lc (Fig. Id). 

First let us make some clarification on the terminology frequently used in the 
literature (hence adopted here too). The contribution to the amplitude of each 
type of correction (" hard- hard" , " hard-soft" and " soft-soft" ) is not uniquely defined 
by itself] only the sum is well defined. The ambiguity essentially arises in both 
approaches from the arbitrariness of the choice of the scalar product (i.e. the 
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weight) with respect to which quantum fluctations A q are taken to be orthogonal to 
the streamline d a A valley , see Eq.(|3^) (valley approach) or to the zero modes zi of 
the second variation of the action around the instanton (in the i?-term approach). 
This problem is often referred to simply as the "constraint dependence" of each 
type of contribution. 

In fact, from the usual insertion of unity in the functional integral, 

1 = J w J d£S{ J d i xw(x)A q (x)f(x)), (40) 

(where / = d a A valley or / = Zi depending on the approach), a fictitious depen- 
dence on the weight is introduced. Such a dependence is expected to disappear 
order by_order in perturbation theory only if all types of contributions are added 
togetherE3'E3. 

Single type of contributions thus in general depends on the choice of weight. 
MuellerE.3 however argued that soft-soft terms are constraint-independent up to or- 
der x 2 . On the other hand, Khlebnikov and TinyakovEj proved (in the i?-term 
approach) that, already at the x 8 / 3 order, the soft-soft contributions to the Holy 
Grail function, F ss (x), are ambiguous: they showed that an appropriate variation 
of the weight induces a variation of the corresponding Levine-YaffeEHI propagator 
of quantum fluctuations giving a modification of F ss to this order. A similar con- 
clusion were obtained in the valley approach by Arnold and MattisE3. Out of this 



observation comes also a suggestive ideacj that if soft-soft corrections exponenti- 
ate but are ambiguous, perhaps the initial state corrections must also exponentiate 
to eliminate the ambiguity. As regards 0(x 8 ^ 3 ) soft-soft term, the constraint de- 
pendence would be removed by hard-soft corrections only, since hard-hard ones are 
known to contribute starting from 0(x 10 ^ 3 ) in the low energy expansion (see below). 

Let us now concentrate on the study of initial-state corrections, considering the 
case of bosonic external particles in the one instanton sector, as often done in the 
literature. The leading order amplitude for a 2 — > n bosonic process can be written 
as 

A 2 ^ n = J dnK( Pl )K( P2 ) J] K(ki)e- S ' (41) 

j=l,n 

where dji stands for integrations over collective coordinates, 



K(p)= lim -m 2 )A mst (p) (42) 

p 2 — >m 2 



are the on-shell residues of the (Minkowski analytically continued) instanton field 
and pi , p2 are the hard momenta. 

The first hard-hard correction comes from connecting together two hard particles 
with a propagator (while other soft lines are kept unchanged), see Fig. 5. It yields 
a correction 

5A 2 ^ n = J dfiD{ PllP2 ) Y[ H(kj)e- 8 ' (43) 

j=l,n 
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where 

D(p,q) = lim (p 2 -m 2 )(q 2 -m 2 )G(p,q). (44) 

G is the Fourier transformed (Minkowski-continued) propagator of the boson fields 
in the instanton background (i.e. the constrained inverse of the second variation 
of the action evaluated at the instanton field, which we call D c ). ft satisfies the 
Levine-Yaffe equation^: 



U c G(x, y) =6(x-y)-J2 fi *j (45) 



arising from constraining the quantum fluctation to be orthogonal to general con- 
straint fi (having non singular overlap matrix f2jj = J dx fj(x) Zi(x) with the zero 
modes of n c , Zj). Clearly, for fi = Zi one recovers the usual BCCL propagators. 

We are thus interested in the behavior of the double residue of the propagator, 
D(p,q), in the kinematical limit of interest pq — > oo, and p 2 — q 2 — to 2 . Explicit 
calculations in the 0(3) two dimensional sigma modelO, and in pure SU(2) gauge 
theory E3 yields 

D(p, q) -c 9 2 p 2 (pq) log(-pq)n(p)K(q) + 0(g 2 (pq) ) (46) 

where R are the residues of the instanton fields, and c is a positive (model depeib 
dent) constant. Applying the Regge-pole tecnique to the operator n c , VoloshirHI 
showed that the constant c is related to the translational zero modes of D c , z M cx 
df,A lst : 

c=(z li ,\x\%)~ 1 - (47) 
Due to the fast growth of D(p,q) with the scalar product (pq) (Eq.p^)), this en- 
hancement overcomes the suppression factor g 2 when (pq) ~ (M\y/g 2 )- The cor- 
rection Eq.(p3|) is no longer suppressed compared to the leading term (an estimate, 
p ~ My/ is used) . 

In a subsequent worktJ Mueller furthermore showed by direct evaluation that, 
for a pure gauge theory, the quantum correction to the two (hard) point Green 
function up to the order are of the form 



N 1 

E ^l(-cg 2 p 2 (pq)log(-pq)) r R(p)R(q) + 0(g 2r (pq)^ 1 )} (48) 



r=l 



where the leading terms of order g 2r (pq) r come from the "squared tree" diagrams of 
Figure 6. The_dominance of the squared trees graph is justified by rapidity ordering 
argumentsO'113. For an alternative (simpler) derivation of this result see the work 
of Li ct al.EZJ. 



If one takes only the leading terms of Eq. ( 48 ) and sums the series up to N = oo 



the hard-hard terms exponentiate into a factor 



e -cg p (p?)log(-pg) ^g^j 
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which gives a contribution of order 0(x 10 ^ 3 ) to the Holy Grail function '. This 
(naive) exponentiation could be taken as a first concrete hint suggesting the possi- 
bility of taking into account the initial corrections semi-classically. 

The first attempt for a semi-classical treatment of initial corrections was done by 
Mattis, McLerran and Yaffe@. Their idea consists of including the initial particles 
into an effective action S e f / by means of the identity: 



J dAA( Pl )A( P2 ) 11 A(k 3 )e lS ^ = J dA\{ Afoje* 5 " 



.(A) 



(50) 



where 

2 

S eff (A) = S(A)-iJ2^gA{p i ), (51) 

i=l 

(the soft fields A(k) can be be treated with the usual ii-term method). The solution 
of the (non covariant) field equation coming from the variation of S e ff, 

5S 



SA(x) 



:E eiPiI Mfe). (52) 



could in principle be used to give a semi-classical approximation that includes the 
leading initial effect for x ~ 1 . 

Another suggestion was made by MuellerfEI. The main idea is to exploit the 
arbitrariness in the choice of the constraints for quantum fluctations, (in the LY 
propagator) to simplify the form of initial state interactions as much as possible. He 
derived the behavior of the generic LY propagator in the kinematical (Minkowskian) 
region of interest ( i.e. {pq) — > oo and on shell) and was able to find an appropri- 
ate (noncovariant) choice of the constraints that eliminates the leading term of 
0(pq log(— pq)) of the LY propagator. This choice consequently eliminates all the 
leading multi-loop hard-hard corrections. As for the hard-soft interactions, Mueller 
demonstrated, by using a diagrammatic analysis that, for a single soft particle the 
interactions with hard particles can be included just by substituting the instantonic 
background by the solution of the (Minkowskian) Y.M. equation with a source 
term: 

Df{W) Gl v (W) = J a v {W) - ±D<?{A) D*(A) (53) 

where D(W), D(A) are covariant derivatives respectively in W or in the instanton 
background, G b fJ ^ l/ (W) is the field strenght in terms of W and £ is the usual gauge 
parameter. The source J contains in a complicated (nonlocal) way the space-time 



^Note that the condition of factorization E/N 2> 1/p (high energies in each final branch of the 
tree graphs) and of exponentiation N 2> g' 2 p 2 (pq) seems to force i<1 

c Also in this case, if the Minkowski propagator had a singular behavior, the loop contributions 
could be important. Approximating the solution by the usual instanton and doing perturbations 
around this (wrong) background, the authors found that a cancellation arises between potentially 
large one-loop terms and that tree contributions reproduce the Mueller's corrections: this suggests 
that in the correct background loops are negligible. 
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integral of the solution W , and has a singularity structure that makes the continu- 
ation to Euclidean time impossible. For more than one soft particles, hard effects 
are to be taken into account with the usual i?-term scheme (all tree graphs should 
be summed) in the new background W. A Euclidean valley resummation does not 
seem to be applicable due to the intrinsically Minkowski singularities. 

It is not clear at present how useful the two approaches mentioned here might 
be in practice, as they (both) involve the solution of an extremely difficult (non- 
covariant) Yang-Mills equation with a source depending in a complicated manner 
on the solution itself f ; they are furthermore incomplete in the sense that they must 
be complemented with some method for resumming soft trees. On the other hand, 
these innovating works perhaps tell us that a semi-classical approach including 
initial corrections is in principle possible. 



2.2. Multiparticle approach 

A more promising approach for estimating the complete Holy GtiaiLfunction (i.e. 
including initial state corrections) is that of Rubakov and TinyakovO'E-X The idea is 
to compute first the inclusive cross section rii — > all for large rii (rii = p-), for which 
a semi-classical method is likely to be applicable straightforwardly (see below). 

2 

Then z/j can be sent to zero afterward in F(x, v{) — logo - >a vj, assuming that 

the limit is smooth, g 2 log 02 ~ Um^^o F{x, Vi). 

Note that while for 02 the exponentiation of quantum (loop) corrections would 
appear rather miraculous, it is quite reasonable that multiparticle inclusive cross 
section can be evaluated semi-classicaly (this will be proved below). 

There are two ingredients for such semi-classical calculations. First, an appro- 
priate Minkowski boundary conditions (taking account of the quantum number, 
energy, etc. of the initial state) must be used. This is to be compared to the stan- 
dard calculation of perturbative S-matrix elements where the vacuum (Feynman) 
boundary condition is used. For this purpose the coherent state representation of 
the S-matrix turns out to be particularly suited. 

Secondly, the fact that one is interested in nonperturbative, classically forbidden 
processes, forces one to consider the time evolution partially in Euclidean direction 
(see below). As a result the classical equation (and its solution) will be defined 
along a complex time contour. 

For completeness, let us first take a few steps back and review briefly the coherent 
state representation for the S matrix (Khlebnikov et alS have used this method 
extensively in formulating the i?-term method.) 

First consider the one-dimensional harmonic oscillator, (generalisation to mani- 
degrees of freedom is straightforward). More details can be found in textbooksGI 
In terms of the vacuum state |0) and creation operator A\ a (non- normalized) 



^One must also check, in the context of baryon number violation, that the solution has nontrivial 
topological number. 
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coherent state can be denned, for every complex number a as 

\a)=e A ' a \0). (54) 

It is easy to derive from commutation relation [A, A^] = 1 the following properties: 

A\a) = a\a) (55) 
(b\a) = e b * a (56) 
e f AlA \a) = \ae f ). (57) 

The wave function of a coherent state in position representation is also useful: 

(q\a) = e -4« 3 -^9 2 +V5^»9 ) ( 58 ) 

(it can be checked that this satisfies Eq. (^5|) when A is expressed in the coordinate 
representation: A — -4= (V^Z + ^j J^))- The action of an operator S in the coherent 
state representation is obviously given by its matrix elements: defining the kernel 
S(b,a) = {b\S\a), 

(b\Sm = [ ^e- aa S(b,a)^(a), (59) 

J 27TI 

where the exponential takes into account the non-normalization of our coherent 
states. Note also that in the integral we must fix a = a* . 

Finally note the following relation involving the matrix element of an operator 
S between fixed number of quanta and its kernel: 

(n\S\m)=d™d£S(b,a)\ a= - b=0 . (60) 

That is all one needs to construct the coherent state_L£presentation of the scat- 
tering matrix (in the Fock's space of asymptotic states) EI Starting from the usual 
interaction representation of S 

S = T \im Jim e iHoT fU(T f ,T l )e' iHoT ', (H = J dk^A^A^ (61) 

and introducing the resolution of unity in terms of fields at Tj, Tf (which allows to 
use the path integral representation for (<j>f\U(Tt , Tj)|^j)), we obtain the following 
expression for the kernel of the S matrix: 



5(6, a) = / dcf> i d ( f> f dcl>e Bi(a ^ i)+Bf(b ^ ) e iS W (62) 
where the boundary terms 

£i(a k ,<k) = -~J rfka k a_ k e _2m k Ti - 1/ dkw k ^(k) ^(-k) 

+ J dky2^^(-k)a k e" iw k Ti (63) 
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+ Jdk y /2^b k ^(k)e t ^ (64) 

come from wave functions of coherent states (see Eq.([58|), and also Eq.(|57|)). 

In the multiparticle approach of Rubakov and TinyakovEl, one actually considers 
another, slightly different kernel: the idea is to project the initial coherent state onto 
the eigenstate of an operator 

= Jdkf k AlA k (65) 

(typically the energy) with a fixed eigenvalue (generalization to cases with more 
than one operator is straightforward). One is then interested in: 

Sn(b,a) = (b\SP n \a). (66) 

Substitution in Eq.(|6^) of the integral representation of the projector 

Pn = J d£e l (°- n K (67) 

and use of property Eq.(|5^), yield 

Sa(b,a) = [d£ [ d&d^f d^e-^ e B ^ a k^ ^ +B i^t) e *sW (68) 



where the function fy. refers to the operator O, Eq.(p5[). 

A quantity of interest, then, is the transition rate for the microcanonical ensam- 
ble (respect to 0). It is obtained by summing the inclusive cross sections over the 
initial states with a definite value of O, 

a(Q)=J2\(f\SP^)\ 2 . (69) 

i,f 

(Actually, there is a slight lack of precision in this expression. One is really 
interested here in a partial sum over B- (and L-) violating final states /, so that the 
unit operator part of the S-matrix actually drops out. Otherwise, Eq.(|69|) would 
have no physics contents! Formally, however, the sum over / is done without any 
restriction (see the next equation). The restriction over selected final states is to 
be taken into account by a judicious choice of the classical backgrounds.) 

By substituting Eq.(|68|) in Eq. (|69|) , one obtains: 

a(Q) = J dcj)d(j)' dadadbdbd£e w (70) 

where 

W = -iQf- J dkb k b k - J dke^k 5 o k o k + Bi(a k ,&) + B f (b k ,(j) f ) + iS{cj)) 
+Bi(a k , #)* + B f tb k , <fi' f )* - iS(<t>')*. (71) 
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Note that a rescaling a^. — > a^e l ^k^, — > aue , the integration over <i(£ + £') 
as well as the substitution £ — £' — > £ have been made. 

The nice feature of cr (fi) is that the functional integral appearing in Eq.(|7C|), 
although very complicated, can be evaluated semi-classically in the limit of small 
g, provided that the external parameters f2 are of the order of p-. Indeed, rescaling 
as 

(cf), </>', a, a, b, b) -> ((/>, <j)',a, a, b, b)/g (72) 

one finds W = tt, where W depends on the rescaled fields, on g 2 fl and on M\\r, but 
not explicitly on g. Thus the usual saddle point tecnique can be applied, yielding a 
result of the form 

ff(fi) = exp{— F(g 2 n,M w )}. (73) 
9 

The saddle point field will be a solution of a classical equation of motion, with 
some new type of boundary conditions (different from those in the usual vacuum 
to vacuum transitions) involving the saddle point values of the complex fields a, b; 
more on this later. 

(The ijresent formalism was actually introduced earlier by Khlebnikov, Tinyakov, 
RubakovlI3, in which the operators O were chosen to be the four momentum. They 
considered the microcanonical ensemble of initial states with a fixed momentum 
(E, 0) , with the idea that the new boundary conditions will be more physically 
relevant in processes of production of many particles, as they keep memory of the 
incoming and outcoming fields. The relevant saddle point configuration is an an- 
alytic continuation to an adequate complex time path of a real Euclidean time 
configuration, the so-called "periodic instanton". a(E) however turns out not to be 
directly related to the cross section oi one is interested inO.) 

The main new idea of Rubakov and TinyakovE! is to consider the cross section 
for the initial microcanonical ensemble with fixed energy E and fixed number of 
incoming particles rij = Vi/g 2 , a(E 1 n i ). 

The advantage of this approach is that for Vi ^ all the initial states corrections 
are automatically taken into account in the semi-classical approximation, since the 
initial particles are represented as fields appearing in the boundary terms. Thus, 
if 02 is recovered in the limit z/; — + 0, it is likely that the initial corrections are 
accounted for: this would reduce the problem of evaluating the Holy Grail function 
effectively to the search of an adequate classical configuration. 

It was also shownLa that, perturbing the functional integral for a(E, ni) around 
the instanton (expansion valid for E <gC E sp ), that the leading hard- hard corrections 
thus inferred do coincide with the naively exponentiated Mueller's series, Eq.(f49|). 

Let's us consider thus the expression Eq.(|7(i|) for a(E,rii) (with the obvious 
generalizations: d£ — * d^d-q, £1 — > -E,n.;, — > wjj, 1), and discuss the boundary 
conditionals that characterize the saddle point field. 

The integration of the final states b, b gives exactly a term 8{<j)f — 4>'f) (by com- 
pleteness of coherent states); then the subsequent integration over <j> (keeping 



R. Guida, K. Konishi and N. Magnoli 21 



into account of all boundary terms g ), gives a 6(<j)f — <pj). Thus, being <p — (p' , <p = <p' 
at Tf = +00 the two classical solutions coincide over all Minkowskian space. 

The integration over a, a is Gaussian and can be done exactly; the integration 
over <pi, (p\ leads to the following boundary conditions at Ti — * —00 (the saddle point 
values for a, a have been substituted already): 

z4(k)+c k 4(k) = e * A k^(k)+^(k)] (74) 
»0i(k) - ^(k) = e^ A k[z^(k)-^(k)] (75) 

where A k = + 77. Subsequent saddle point evaluation of the integral over £ and 
77, will then fix the parameters of the solution in terms of E,rii. 

By time traslation invariance, the variable £ associated with ejqem y can always 
be taken as purely imaginary: £ = i£o; perturbative calculationsL3'Lj suggest that 
also the saddle point valufLpf 77 is purely imaginary; so we can choose 77 = ir/Q. It 
was furthermore observedE3 : L3 that the parameter £0 can be removed from boundary 
by choosing suitable contours in complex time plane: one 
choice could be the one shown in Fig. 7 (with £ = BB'). 

If we assume the solution to be free at r = Ret — > -co, and write in this limit 



0(k) = / k e ^k T + f_ k e l "k T teAB (76) 

<£'(k) = 5k ^k r +5_ k e^ teA'B', (77) 

substituting in Eqs.( |74| - |75| ) we obtain simplified boundary conditions: 

fk = e ~ V °9k; Ik = eV °~9k (78) 



which relates asymptotic positive and negative frequency components on the two 
branches of the path, AB and A'B'. 

As anticipated before, the necessity of using the complex-time deformed path 
comes out from the double requirement of having free particle incoming (outcoming) 
at Minkowskian time Rei — » — (+)oo and of describing a semi-classical, classically 
forbidden transition, that requires a part of the evolution to proceed in the Euclidean 
time. 

The idea of the complex-time deformed paths is not really new: it is known 
(and was invented first) in quantum mechanics, in the W.K.B. anproximation for 
the fixed energy Green's functions in the case of barrier penetrationEi In that case, 
the relevant classical solution is composed of Minkowski solutions for the initial and 
final parts of evolution and of an intermediate, Euclidean solution describing the 
tunnelling. 

SNote that also a term 

5 / <* k M-k-^'-k^f 

in the exponential coming from integration by parts of the kinetic term of the fields' action, are 
essential to the integration over initial and final fields. 



conditions Eqs. ([74-75 
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Let us now have a closer look at the required properties of the classical field 
we are looking for. If the solution to field equations satisfying boundary conditions 
Eq.([?8|) is unique, then it must obey the relation 4>(t*) = 4>*(t) (coming from reality 
of field equations); that implies the reality of the field on the positive Minkowski 
time. It follows from Eqs.(|7§-^) that 

fk = 9jf\ 5k = /£• ( 79 ) 

Also, it can be seen that the saddle point values of a, a are complex conjugate of 
each other, which means that the sum over the initial states is dominated by a 
single coherent state. 

Perturbative calculationsOCj support the conjecture that on the Euclidean time 
piece the solution might be real: this (together with the reality on the positive real 
time axis) would imply the presence of a turning point = for t = OlJ. 

All these hypotheses, put together, imply that the required saddle point config- 
uration, determining the microcanonical cross section a(E,rii) would be a solution 
of classical equations on the complex time path of Fig. 7 with boundary conditions: 

0|t=o = 0; 

0(k) = f k e-^ T + e^r k e^\ (80) 

where r = Re t — > — oo. 

Unfortunately, the "microcanonical" saddle point is not known at present even 
numerically at energies of order_E sp . A perturbative expression for this solution 
(valid at low energies) is known in the context of two dimensional Abelian Higgs 
model: it looks like a chain of alternated instantons and anti-instantons. This 
study put into evidence the problem of the presence of singularities of analytically 
continued classical solutions in the complex time, that the chosen path must not 
touch. See the discussion below, Eq.(p^). 

A slight variation of the approach of Rubakov and Tinyakov, is to compute the 
cross section for a given coherent state of fixed energy 

cr(E,{a k }) = J d(j)d<j)' dbdbd^d£_' e w (81) 

with 

W = -iEfc-g)- /dk6 k 6 k + J B i (a k e ia 'k^ i ) + J B / (6 k ,0 / )+ 1 t5(0) 

+B,(a k e ia, k«' ) ^)* +£^,4)* -iStf)*, (82) 

where now a, a are the arbitrary but fixed complex numbers. Integration over final 
fields goes as before. What changes is the integration over 0j,<^ that by saddle 
point evaluation give the initial boundary conditions: 

*0i(k)+^(k) = ^/2^a k e" ia, k (T *- s) 

ito-cVKk) = -y/2^a_ k eV*- e) . (83) 
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Clearly, if one would have known the solution of the boundary condition problem for 
generic values of a, a one could estimate semi-classically Eq.(|Sl|) and recover directly 
the two particle cross section by differentiating twice o~(E, {a^}) (see Eq. (|(io|) ) . 

Unfortunately, the resolution of the general boundary condition problem is a 
hard task. Some recent work EBB pursue a more modest goal of calculating semi- 
classically o~(E, {a^J) w ith a i o fixed a posteriori to satisfy the boundary conditions 
Eq.(^3|), starting from some given classical solution. These works have many com- 
mon features. They study exact solutions of field equations in conformally invariant 
models 11 , obtained by some symmetry requirement thaLsimplifics the problem to an 
equivalent one dimensional one. In particular, the firstlla studies the massless two <iL 
mensional 0(3) sigma model, searching for 0(2) symmetric solutions; the secondlZa 
considers the 50(4) conformally invariant solutions of Minkowskian SUCZ) pure 
Yang Mills theory, the so-called Liischer-Schechter solutionsSEa; the thirdE3 stud- 
ies an 0(4) invariant solution of a massless, four dimensional 4> 4 theory with the 
positive ("right") sign of the coupling constant. 

All these solutions are continued to a complex path as in Fig. 7, and possess a 
turning point at t = 0. It turns out that the position of the path relative to the 
singularities of the solution determines its very nature. As a significant example 
note that Yang Mills complex time solution is such that the Minkowskian action S 
and the topological charge Q evaluated along the path, obey the relation 

I^ImS = Q = N, (84) 

where N is the number of singularities between the path and Minkowski time axis 
(relation that is identical to the usual multi-instanton). The path considered in his 
particular case corresponds to the one instanton sector. 

Evaluating the average number of initial (from the derived a, a) and final par- 
ticles (from the saddle value of b,b), it turns out that these solutions allow to 
compute the probability of a process with parametrically less particles in the initial 
state than the final one (for example^ Vf — vj^ 8 for small Vi). But in the limit 
V{ — > the cross section is exponentially dumped by the full t'Hooft suppression 
factor exp{— \6ir 2 /g 2 }. These solutions however do not maximize the transition 
probability in the one instanton sector at a given energy, (because the resulting a, a 
have not opposite phase, see above). Thus these solution neither are directly related 
to the 2 — * all processes, nor can be used to establish a rigorous upper bound to 
that cross section. Nevertheless they are a useful benchmark to understand the 
property of more physically relevant solutions (such as the "microcanonical" one 
considered before). 

An important, nontrivial question, is if general complex gauge fields configu- 
rations defined on a complex contour, satisfy a generalized Atiyah-Singer index 



^They are all massless models that are supposed to describe the high-energy dynamics of the 
massive ones. Clearly with this approach it is impossible to recover the low-energy behavior of 
massive models. 
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thcorcmli2jiiij of the form 



m-n r = / TtFF = Q (85) 

J path 

where ni (n r ) are the left (right) handed zero modes of the Dirac operator. Such 
a theorem would imply that if we had a classical solution with complex topological 
charge Q — ±1, we would be assured the presence of the fermion zero mode which 
is essential for the nonvanishing of the anomalous cross section. 

Also, one may ask if there is a level of the Dirac Hamiltonian crossing zero in 
time evolution (physically interpretable as creation or annihilation of a hole or a 
particle depending on the direction, see Section 5), in one-to-one correspondence 
with the zero modes. A step in the direction of answering these questions, was taken 
in a work by Rubakov and Semikoz LD, where the authors verified these conjectures 
in the case of the two dimensional Abelian Higgs model. 

Let us end this section by a brief mention to a recent work of Muelleio. Instead 
of the " microcanonical" cross section <r(E, vi) of Rubakov and Tinyakov, the author 
considers a slight variation, with the initial state consisting of the (equally weighted) 
sum of all pure coherent states with energy E and with n\ = v\j ' g 2 right moving and 
n-2 = v 2 /g 2 left moving W's (he considers only the pure gauge sector). The processes 
with such modified initial states seem to be considerably closer kinematically to 
those with the two particle initial states, than the simple microcanonical ones. 
Subsequently the functional integral is expanded perturbatively around the single 
instanton. An accurate study of the hard- hard corrections, (up to x 10 ^ 3 in F), and of 
hard-soft (up to a; 16 / 3 ) shows that in the limit v\ — v 2 = v i —> these contributions 
to the exponent F(x,Vi,V2) join smoothly with the result for two particle initial 
state. This is quite nontrivial, and holds only after compensations of divergent \/v 
terms. 

It would thus appear that a semi-classical (perhaps numerical) study of processes 
of the form n\ + n 2 — > all provides a good way to estimate nonperturbatively the 
Holy Grail function for the desired 2 — > all transition. How realistic this possibility 
is in practice, is however unknown. 



3. Quantum mechanical analogue problems 

In order to get an intuitive, physical understanding of our four dimensional prob- 
lem, several authors studied toy models in the context of one-dimensional quantum 
mechanicsE3 _ EJ. The results of these studies will be reviewed briefly in this section, 
following mainly the treatment of Diakonov and PetrovO. 

Consider a quantum mechanical double well coupled to a weak, rapidly oscillat- 
ing field. The Lagrangian is given by 

L = 7 Q ( il )2 - l {q2 - y2)2 - Fq{elEt + e ~ lEt) ) ' (86) 

where F is a small coupling constant. One could also introduce fermioiiic excitations 
(instead of the external field), without however changing the resultsEJ essentially. 
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Let us ask now the following question: what is the probability of the transition 
from the ground state in the left well to a highly excited state in the right one? To 
lowest order in the external field, the transition probability from the ground state 
to the iV-th state is given by the well known Fermi's Golden Rule: 

W = 2it5{E + E - E N )\(N\q\0)\ 2 F 2 (87) 

where 

(N\q\0) = [ dqV* N (q)qy (q). 



The problem is then reduced to the calculation of the matrix element of the position 
operator betwen two states with a large difference in energy. 

Two approaches have been taken to estimate this matrix element: one uses the 
semi-classical (W.K.B.) method (which will be followed below) and the other (taken 
by Bachas) a more formal one. 

As is well known, the nonperturbed energy eigenstates of the double well poten- 
tial are symmetric or antisymmetric under parity. For a level lying much below the 
height of the central barrier, one can define states localized in the left or right well: 

\Nl) = -j=(\N.) + \N a )) (89) 

Wr) = ±(\N.)-\N a )). (90) 

It is clear however that when the energy is higher than the top of the barrier these 
states are no longer localized in the left or right well. In the Weinberg Salam the- 
ory in four dimensions, the states corresponding to different Chen-Simons numbers 
are always well defined (for instance they can be distinguished by the number of 
lefthanded fermions present). Thus the analogy between the real problem and the 
quantum mechanical toy model is not self evident for energies above the barrier. 
But the matrix clement is exponentially small in such a situation (see below) so 
that the precise definition of the final state is probably not crucial. 
The matrix element is given, to exponential accuracy, byE3 



where 



(N R \q\0 L ) = e (-*+*-*), (91) 



S 1 = -L / dq(q 2 -v 2 ), (92) 



V J- 



1 flo 

S 2 EE — / dq^{q*-v 2 ) 2 - 
%9 J-oo 



(93) 



S, ^ I <h lvU f-ri)l-<.. (<I4) 



<IF 



-,2 



where e ee 8Eg . The terms in the exponential are shortened action with energy 
or E, corresponding to the (complex time) trajectory shown in Fig. 8. To our 
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purpose the only important contributions come from the parts of the trajectory 
corresponding to Euclidean time evolution (AB, CD and EF in Fig. 8), because 
the paths in Minkowski time yield only a phase factor. Note that when E — > the 
contributions from paths AB and CD tend to cancel each other while the path EF 
gives the usual tunnelling factor at zero energy. 

As the energy increases, it is clear that 5*3 diminishes but Si — 5*2 grows with 
energy. The problem is which tendency dominatesE3. 

The answer has been found for very high and very low energies. Furthermore, an 
argument has been given for the "half-suppression" i.e, the transition probability 
at its maximum (which occurs at E — v /8g 2 ) being the square root of that at 
E = oB@ 

The transition probability W is proportional to the square of the matrix element 
(iV|q<|0). It can be calculated without an explicit knowledge of the trajectory at very 
high or very low energies. Indeed, first consider the quantity 

2 d\ogW 2 dSi 2 dS 2 2 dS 3 / os \ 

9 ^^-' 9 ^T + 9 ^T +9 ^T- (95) 

The right hand side can be simplified: if e < 1 it is 

X _ i _£!*> (96) 



while if e > 1 it is equal to 



= ^E, (97) 



de 2^/V^ 



where K(x) is the complete elliptic integral of the first kind0, k = \J \ + ^ and 

After integration in e one finds the following results: 

5 2 logW^-^ + ^(log^ + l) + 0(e 2 ) (98) 



for e << 1; 

6 v 4' 2 



/i 0# = _i ( li ) + ( £ V^ (99) 



for e >> 1, where B(x,y) is the Euler Beta function@. 

These results show that the transition probability is exponentially suppressed 
both at low energies (where W grows with energy) and at high energies (where W 
is exponentially damped with energy). 

To demonstrate the square-root suppression of the transition probability at its 
maximum, one compares the (imaginary) time along the paths CD and EF. Ac- 
cording to the classical mechanics 

fl a ^=T 7 (e), (100) 
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where Tj(e) is the time along the 7-th path with energy E. The time interval T(p/) 
of the piece of the contour 7 is thus given by 

™-lvw=< (101) 

Let us consider the difference of time intervals along the contours 71 and 72 (Fig. 
9), in the complex q plane. As long as the potential has no singularities for finite 



values of q, the integrand of Eq.(lOl) has only cuts at the classical turning points. 
The cut being the square-root type, the value of the integrand at both sides of the 
cuts differs only in sign. Using Cauchy's Theorem one gets then 

T 1 +T 2 = 2T 3 (102) 

or 

_ds 1 = 2 ds 1 _ 2 ds 1 

de de de 

and after integration in e, this yields 

S 3 (e) = S 3 (0)-2S 1 (e) + 2S 2 (e). (104) 

At e = 1, 53 vanishes, thus 

S 1 (1)-S 2 (1) = ±S 3 (0). (105) 

This result connects the overlap integral (the left hand side of the equation) to 
tunnelling factor at zero energy, as anticipated. This is a very general resultO 
independent of the explicit form of the double well used above^ 

Essentially the same result has been obtained by BachasEa, who was able to 
give a rigorous proof of exponential suppression of the induced high-frequency tran- 
sition amplitudes in an anharmonic potential (one or two wells), based on exact 
recursion relations between matrix elements of powers of the position operator. 
This demonstration settles the issue for the quantum mechanical problem, but, in 
spite of optimistic expectations of the author, no generalization to field theory is 
known, up to now. Finally we find it interesting that the high energy behaviour 
of the "Holy Grail function" (for the case of the double well), Eq.(p9|), found by 
the WKB method, is consistent with the exact bound found by Bachas, apart from 
a logarithmic term which is probably beyond the precision of the semi-classical 
approximation. 

Landau's resulL-Eq. (^Tj) , has been earlier generalized to field theory by Iordan- 
skii and Pitaevski£3. The asymptotic behavior of the imaginary part of the Fourier 
transform of a two point Green function is given by: 

Im G(k, u) ~ e 2 (- AS '(°<°)+ AS "^ k » (106) 

The two terms in the exponent correspond to the value of the action along the 
trajectories at energy E = and E = to. This formula is valid if u> and k are large 
quantities, and is similar to the quantum mechanical one. 
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The method of Iordanskii and Pitaevskii was recently applied to non Abelian 
Gauge theoriesO. 

By analogy to the quantum mechanical example one must find singular solutions 
to the classical Yang-Mills-Higgs equations of motion at energy and E. Diakonov 
and Petrov look for an 0(3) symmetric solution in the pure gauge theory. At 
energy these are the instanton and anti-instanton with the size p 2 changed to — p 2 , 
while it is more difficult to find a solution if the energy is non zercO. In the limits 
of very high or very low energies an approximate solution can however be obtained 
analytically. In the low energy regime, the two contributions to the reduced action 
corresponding to the path ABCD and EF in fig. 8 sum up and give the usual leading 
result Eq.(|J). 

In the high energy limit, the cross section for isotropic multiparticle production 
is found to decrease exponentially with the energy as 

a(E) ~ e -const.( ap Efro /a ^ (1Q?) 

(valid for fixed p such that apE ^> 1.) This result can be interpreted^ in terms of an 
effective strong repulsive interaction between the " instanton" and " anti-instanton" 
at small separation. Such a repulsion prevents the instanton pair from collapsing, in 
contrast to what happens in the case of the valley (Section 5). (See Klinkhamertl 
for a related discussion.) 

The result should be correct also in the electroweak model, because both at very 
high energies and low energies one can neglect the effects of Higgs boson, while at 
the sphaleron energy one should consider the complete Yang-Mills-Higgs coupled 
equations. 

A naive extrapolation of the results at intermediate energies shows that the 
maximum of the cross section is achieved near the sphaleron energy where it is close 
to the square root of its value at zero energy. One cannot draw any conclusive answer 
because it is not clear whether the particular classical solution used maximizes the 
multiparticle production cross section. 

Finally the work of Cornwall and TiktopouloSEI should be mentioned. These 
authors use a functional Schrodinger equation approach to study SU{2) gauge the- 
ory, exploiting the analogy with the quantum mechanical double- well problem. First 
it is argued that the matrix element, 

(NE\m\0) = J #^{0WoM (108) 

(where the final state with energy E consists of N particles, <j> is the field operator), 
is bounded from above simply by 

M4>ne}- (109) 



4>ne is the appropriate field configuration with quantum numbers E 7 N. Secondly, 
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the vacuum wave functional for the SU(2) gauge theory is knownEil (</> = A): 

2vr 

exp - 

J— — oo 



^o{A}^AA ex P 



-\W(A)-J\ 



(110) 



where W(A) = fifes is the Chern Simons functional (see Eq.( |17§| )). By approx- 
imating the final state configuration A by the sphaleron for which W(A) = 1/2, 
and by keeping the J — 1 term only, these authors get the half (or the square root) 
suppression of A(B + L) ^ cross sections. 

4. Unitarity Bounds 

Various arguments based on unitarity, which all lead to the "half-suppression" 
of the A(B + L) cross sections, are reviewed in this section. 

4.1. Multi-instanton Militarization, "half-suppression" , "premature uni- 
tarization" and all that 

The result of the original calculation of RingwaldB and EspinosaB violates 
unitarity if extrapolated to high energies, for reasons similar to those for which 
the Fermi theory of weak interactions violated unitarity The effective interaction 
vertex is given by a local, non-renormalizable formt3li2!J't3: 



L 



= I dx 



dp 



du d(p)e 



-2-rr 2 p 2 4>(x)4>(x) _ 



(111) 



where U = U^a^, U = U^a^ are the matrices representing a global orientation of the 
instanton in the isospin space, and the instanton density d(p) contains exp(— S mst ) 
and the functional determinant. 



The use of Eq. ( 111 ) in the lowest order to compute the S matrix elements will 
necessarily lead to violation of unitarity. Contributions which are higher orders in 



Eq. ( 11_1 ) - multi -ins t ant on contributions - have to be taken into account in order to 



aptH 



restore unitar 

ZakharovE£3 (also Aoyama and Kikuchll2il, Veneziandl23) pointed out that such 
multi-instanton corrections can change qualitatively the result found in the original 
one instanton calculation. A rough estimate of the multi-instanton contribution to 
the 2 — > N amplitude, iterated in the s-channel (see Fig. 10) and summed over the 
number of instantons, would give an answer, 



A 



(j A tree \2k 



E°° A tree 
fe=0 A 2^I 

~A^ n /{1 + (A%Zn) 2 }, 



(112) 



where a simple factorized form of k- instanton contribution is assumed, the super- 
script "tree" indicates the single instanton (or anti-instanton) approximation, and 
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N stands for the dominant, multiparticle intermediate states (it includes symboli- 
cally also the summation over them). Although such a simple, factorized structure 



of Eq.(112) is by no means obvious, a more careful estimation of multi-instanton 



contribution will be given in the next subsection which will reproduce essentially 



all the results following from Eq.(112) 



The single instanton amplitudes can be taken, for the purpose of resummation, 
to be of order, 

A l 2 Z e N -e- 27r / a -K 1 / 2 , (113) 

and 

A%™ N ~ e- 2lx l a ■ K, (114) 

where K is the factor due to the sum over soft gauge bosons produced by instantons. 
At low energies K is known to grow exponentially as, 

1,2. f E -,4/3 

K~ e + ; c=O(l)>0. (115) 



The difference between Eq.(113) and Eq.(114) reflects the fact that in the latter 
there is a sum over both the initial and final (multiparticle) states while in the 
former the summation is only over the final states. 



It follows from Eq.(112) that the unitary amplitude is limited by 



RK 1 / 2 

A 2 ^N = 1 + B2R2 < const. B 1 ' 2 , B ee e" 2 ^, (116) 

hence 

cta(s+l)#o ~ |^2^jv| 2 < exp , (117) 

the well-known "half-suppression" result. The terminology is due to the fact that 
at the tree level oa(b+l)^o ~ ex P — "S" 



Note that this bound (Eq.(116), Eq.( 117 )) is independent of the way K depends 



on the energy, i.e., does not depend on the low-energy approximation, Eq. flllq ). It 
is thus probably of a broader validity than it might first appear to be. On the other 
hand, there is no guarantee that the upper bound is actually reached: the correct 
statement is that 0a(b+l)^o is suppressed at least by the half (or more properly, 
the square root) of the standard 't Hooft factor. 



If one does use the leading semi-classical result, Eq.(115), a somewhat stronger 
result follows. Namely, the full amplitude would be dominated by the multiinstanton 
terms in Eq. ( |112| ) , as soon as the energy reaches a critical value such that A t ^ e f tN ~ 1 
(it turns out to correspond to the value E max ~ 0.95E sp , in the leading semi-classical 
approximation: see Eq.( |139| )). But at this energy the tree amplitude A\ r ^ N is still 
exponentially small! 

It would appear then that the multi-instanton corrections overwhelm the single 
instanton contribution at an energy such that the latter is still exponentially small, 
invalidating any argument based on single instanton calculation extrapolated to 
higher energies This scenario was termed "premature unitarization" by Maggiore 
and ShifmanEil 
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At preserraii3 t!3, there is no rigorous proof of premature unitarization, as for- 



mulated above. The estimate Eq.(115), based on the leading semi-classical approx- 
imation, may receive important modifications when higher order corrections are 
taken into account. However, the alternative possibility - that the multi-instanton 
contributions never overcome the single instanton term - means that the K factor 



appearing in Eq.(112), Eq.(116) is never able to compensate the original t'Hooft 
suppression factor. In other words, the single instanton calculation may be saved, 
but it must stay exponentially small, after all! 

This discussion illustrates the essence of the arguments based on unitarity. In- 
dependent of the details (for instance, whether or not premature unitarization holds 
true), multi-instanton unitarization implies the exponential (at least "half-") sup- 
pression of cr A(s+L) ^ - 

A similar conclusion, in fact, follows also from a purely Miukowslri picture of 
our processes: the sphaleron resonance formation and its decayu'ESlIia. 

Consider the Breit-Wigner formula for the 2 — > N amplitude 

r l/2 r l/2 

{E - E sp ) + iT tot ' V ' 

where T to t — is the total decay width, T 2 is the partial width to the two particle 



state. Eq.(118) satisfies unitarity, Im^2^2 = J2n l^2^w| 2 . Now, the sphaleron is 



a classical field with spatial dimension ~ 1/gv, hence the decay particles will have 
momenta of the order of ~ gv . It follows that the average number of the secondaries 
is N ~ E sp /gv ~ 1/q (recall E sp ~ v/g.) If one assumes a Poisson distribution for 
the number of the secondaries - as is appropriate for a coherent state - then one 
ends up with the estimate 

\a 2 ^ n \ 2 < ryr tot ~ e ~ R ~ e - const >. (119) 

Thus the cross section for the sphaleron formation is small because the latter is 
coupled weakly to the two-particle initial state. 

The fact that we find from the sphaleron picture an estimate similar to the one 
obtained in the multi-instanton picture, is not really surprising. The sphaleron is 
an unstable state staying on the top of the barrier separating two adjacent vacua 
of the SU(2) gauge theoryEl, decaying with equal probabilities to A(B + L) ^ 
and A(B + L) = final states. In the context of Euclidean, instanton description 
of A(B + L) 7^ transition, the sphaleron should correspond to multi (infinite) 
instanton configurations. 

Finally, it is possible to understand both types of arguments based only on a 
very general feature of unitarity and the particular aspect of involved dynamics - 
dominantly multi-particle production. 

Considered the s-channel unitarity equation for the forward elastic amplitude 



lmA el = \ A i-> 



Nl 2 



N 



32 Fermions, anomaly and unitarity 



= L4 ei | 2 +^LW| 2 

(X (T e l+(ri n el, (120) 

where we assumed a simple form of unitarity appropriate for the S-wave0, that is 
adequate for processes described by a single instanton or sphaleron, and neglected 
all complications arising from spin. As is well known, it leads to the upper limit for 
an S-wave amplitude. 

\Aa\<l. (121) 

The question is whether such a limit can be actually saturated in the case of 
instanton-induced multi W boson productions. 

It should be recalled that in the case of high energy hadron hadron scattering 
(due to strong interactions), more and more partial waves come into play as the 
energy increase s, an d this leads to a much less stringent unitarity limit, the so-called 
Froissart boundHB for the total cross section, 

a tot oc log 2 s. (122) 

(which amounts to |^4 e ;(s,0)| < slog 2 s. ) 



Whether the limit (121) can be reached, depends on the dynamics, even if the 



reaction proceeds indeed through the S-wave only. Let us define the inelasticity r, 
r = a mel /a el = ^ \A 2 ^ N \ 2 /\A e i\ 2 (123) 

N^2 



Eq.(120) and Eq.(123) together imply a more stringent bound, 



0a(b+i/)^o < atot ^ -lraA e i < -— — . (124) 
v ' s s 1 + r 

We see that whatever mechanism (in pure S-wave) gives r 3> 1 would lead to 
VA(B+L)jto much smaller than the geometrical value (~ i). 

In a more careful treatment one should distinguish the baryon number violating 
from baryon number conserving intermediate states in Eq.( |l20| ). The argument 
essentially goes through without modification, though, the reason being that the 
particular mechanism under study - multi-instanton sum - will give the same con- 
tribution to A(B + L) 7^ and A(B + L) = processes, at its maximum (see the 
next subsection). 

The s-channel iteration of instanton anti-instanton chain is a crude approxima- 



tion which ensures the s-channel unitarity, Eq.(12C). A four point amplitude must 
satisfy the t-channel (as well as the u-channel) unitarityc£3 also. That requires at 
least the instanton chain iterated also in the t- or u-channel. Such contributions 
can modify the S-wave dominance of the amplitude: cannot the electroweak high 
energy scattering then become of multi-peripheral type at high energies, with a 
much larger cross section, and become similar to the "soft" hadronic processes? 
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Such a possibility cannot in principle be excluded. Nevertheless, there is a 
difference here as compared to the multi-instanton contribution iterated in the direct 
(s-) channel. The large energy (large s) carried by the initial particles is converted 
into the energy of multiple gauge or Higgs bosons. The effect of the final state 
summation is to enhance drammatically the leading instanton amplitudec2l'E!l. The 
effect is so large as to lead to the violation of s-channel unitarity: one is forced to take 
into account higher order (multi-instanton) contributions, iterated in the s-channcl. 
On the contrary, no reasons are known to suspect that multi-instantons iterated 
in the t-channcl (for instance) are of any particular importance. Multi-instanton 
chains iterated in the t or u channel presumably remain negligibly small. 



4.2. Resummation of multi-instanton contributions 

These arguments above, leading to the " half-suppresion" , based on a very rough 
estimate, Eq. (|ll2| ), will be corroborated below by actually resumming the multi- 
instanton £antributions. The treatment closely follows that by Musso and one of 
the autora! 



The starting point will be the one-instanton i?-term calculation of the A(B + 
L) =/= cross section, or the imaginary part of the forward elastic amplitude, re- 
viewed in Section (1.3). To generalize the formula to multi-instanton cases one 
needs the knowledge of the amplitude itself rather than its imaginary part, since 
the former (more precisely the S-matrix element) can be iterated in the sense of 
Feynman diagrams. 

To compute the amplitude, we recall that if the instantons and anti-instantons 
are far apart, the single instanton S-matrix elements can be treated as an effective 
Lagrangian, L inst . The result Eq.(p8[), Eq.(p9|), Eq.(|30|) just corresponds to the 
lowest contribution in Li nst to the amplitude of W-boson production, squared and 
summed over the final states. Thus to get the amplitude, one replaces the "cut 
propagator" 

" k 2tt / d 4 k9(k°)5(k 4 ) (125) 



(126) 



appearing in Eq.(|2^) by the uncut propagator 

d 4 ki 



- ic 



One finds thus, 

<p,-p\A\p,-p>=ie- 16 ' r2 ^ 2 J dtdxdpdp'C(p,tp',Oe W ' ■ (127) 



The exponent W is given by 



d 4 ki 
fc2" 



-e lkx (3k 2 +k 2 ). (128) 



le 



W> = -lEt - ir 2 v 2 (p 2 + p' 2 ) + ^p 2 p' 2 J 
Now, since 

r j4i„a 4 2 

_yk X = _ s ( 129 ) 



/ 



k 2 + 
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the integral in k in W' is equal to 

327r 2 (3a;2+x 2 ) 



(x 2 — ie) 3 



(130) 



Computing the x integration by the saddle point method, one finds at the saddle 
point (x = 0), 

q fi 2 2 a 

W' = -iEt-n 2 v 2 (p 2 + p' 2 ) + mV fJ =W: (131) 

g t 

namely, in the leading order we find the same exponent for the forward amplitude 
as for its imaginary part. This is consistent with the fact that the leading i — a 
contribution to the forward amplitude is purely imaginary. Computing the resulting 
integrals over t and p, p' by the saddle point method as before , one finds thus 



< p, -p\A\p, p >= i exp( — + 3( Q ^ 2n|4 ) 1 ) , (132) 



167T 2 | g( 3E 4 ,1/3 

g 2 ^8ir 2 g 2 v 

in accordance with Eq. ( |l8[ ) . 

The contribution from n pairs of i — a (Fig. 10) to the forward amplitude can 
be written down by generalizing Eq.(127), Eq.(12J). If we neglect various pre- 
cxponcntial factors, it reads, 

... / 1] duY[d Pl expW. (133) 
J i=i i=i 

2n-l 2n Q „ 2 2n— 1 

W = - lE J2 U - -V E^ + ^r£ PiPi+^t' ( 134 ) 

i=l i=l " i=l 



In arriving at Eq.(134) the integrations over the relative instanton orientations have 
been done by the saddle point approximation which effectively set all instantons 
and anti-instantons alligned in the SU(2) space. The integration over the relative 
instanton (space) positions yields 

xi = x 2 = ... = x 2 „_i = 0. (135) 

The remaining integrations over the relative instanton time coordinates ti, and their 
sizes pi, can also be done easily. See Appendix A. 

The resulting 2n instanton contribution to the forward // — > f f amplitude is 
now (in terms of the dimensionless variable x = E/E sp ) 

A (2n) ^ . 2 „-l i/ 2r 1 -l e ^^((l-l-0817^/ 3 )r l +0.7672^/ 3 )^ 

where the real factor H contains the power contribution from the fermion zero 
modes in the alternative i — a " bonds" as well as the preexponential factors arising 
from the gaussian integrations. 
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Note that the phase factor i 2 ™ -1 in A^ 2n ^ originates from the continuation from 
2n i — a centers in the Euclidean space to Minkowski spacetime positions. (One 
missing i is due to the standard definition of the amplitude.) All other Gaussian 
integrations over U,s and pi,s are real. As a result the sign alternates in the sum 
over n (see Eq.(|112|)), and gives 



oo TT ~ 16 / 2 [l-1.0817a; 4/3 ] - 
A ^y A (2n) ^ l ^[0.767^] 

This concludes the calculation of the forward elastic amplitude in the leading 
semi-classical n-instanton approximation, resummed over n. 



There are several noteworthy features in Eq.(137). 
(i) As expected, the resummed amplitude behaves in a way qualitatively different 
from the single instanton contribution (n = 1 term only). The resummed amplitude, 
which is unitary, turns out to be exponentially suppressed at all energies. It reaches 
the maximum 

\A\ max ~exp- — (1.333) (138) 
a 

at 

E max ~ 0.95£ sp . (139) 
where 1 — 1.0817a; 4 / 3 ~ 0. At higher energies it is exponentially damped. The 



suppression Eq.(139) is somewhat stronger than the naive "half suppression" factor 
exp(— 2n/a), in agreement with our general conclusion. 

(ii) That the resummed amplitude A = Y^=i A^ 2n ^ satisfies the s-channel unitarity, 
(A 2 ^ 2 - A* 2 ^ 2 )/2i = E A 2^n x Al^ff, (140) 

N 

can be seen as follows. Consider A 2 ^ 2 and A 2 _+n, both computed in the leading 2n- 
instanton approximation, summed over n. Compare the two sides of the unitarity 
equation, namely lmA 2 ^ 2 and ^2 N \ A 2 ^n\ 2 . The latter can be written as 



E E 4l N (A m J N r - E E E J&P^r- m) 

N l,m=l N n=l 1 = 1 



Recall that the sum over the intermediate states N (corresponding to cut prop- 
agators) gives rise to the same exponent as the sum over virtual intermediate states 
(with uncut propagators), as noted before (after Eq.(|131[)). Taking into account 



the appropriate phase factor the (n, I) term of Eq.(141) is 



J2^- N l) (4l N r = i 2n - l -\-iT'\A^ 2 \. (W2) 

N 
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These terms are seen to be in one-to-one correspondence with various cuts of multi- 
instanton chain for A^^- Indeed, the identity, 

2n-l 

i Jh *- , - 1 (_i) 1 - 1 = = (i 2 "- 1 - (~i) 2n - 1 )/2i (143) 

1=1 



shows that Eq.( |l40|) is satisfied in this approximation. 

(iii) Note that both types of cuts corresponding to A(B + L) ^ and A(B + 
L) = intermediate states appear in the unitarity equation. However, it is not 
difficult to separate the two types of contributions in Eq.(137), once one identifies 
the contributions of various cuts, Eq.( |l43[ ). One finds 



Im A = Disc A{B+L) ^ A + Disc A{B+L)=0 A. (144) 

-[l-1.0817x 4/3 ] 



Disc A(B+L) ^ A = _ - e « , (145) 

(l + ^2 e -?- [1 - 1 - 0817a/I ) 2 

^[1-1. 0817a; 4 / 3 ] 



H e n' 2 1 ' ' _ i6x£[ 7672x 4/3 l 

Disc A ( B +n=o A = —= e s- , (146) 

(+)0 (l + ^ e ^[l-1.0817xV3 ])2 

and 

o-A(s+i)^o oc (-)£>isc A(B+L)#0 A. (147) 

Clearly all the point made above in (i) applies to &a(b+l)^o which is the quantity 
one is really interested in. 

Finally let us note that at the energy Eq.( J13S| ) at which <Ja(b+l)^o takes the 
maximum value (and where 1 — 1.0817a; 4 / 3 ~ 0), the baryon number violating and 
conserving cross sections are the same: 

CA(B+L)^0 — °A(B+L)=0 (148) 

to exponential accuracy, in accordance with the sphaleron resonance formation pic- 
ture. 

(iv) Instead of 2 — > TV amplitude considered above, one can study N — > N', with the 
coherent states method of Khlcbnikov et alEa. Consider in particular the forward 
amplitude, 

// + N gauge bosons — ► // + N gauge bosons, (149) 

in the multi-instanton approximation, summed over N. The R-term technique for 
doing the sum over the number of particles in the intermediate states used above 
can be easily extended to the initial state (see also Section 2.2). 
One finds that 

\ A ff+N^ff+N\N=N mam ^ | ^2 A ff+N^ff+N\ (150) 

all N 
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for some N max of order of 0((^)(-gP-) 4 / 3 ) by using the well known relation between 
a coherent state and a state with a given number of particles. Also, the right hand 
side of this equation can be computed just as in the case of the forward 2 — ► 2 



amplitude, Eq.(137). The crucial difference is however that now no "end point 



108 



effect" is present and as a result one gets 

^-ff+N gauge bosons — >ff+N gauge bosons | N— N; 



r^_[i_i.0817a; 4/3 ] 
Hen 21 

1 + H 2 e ^[l-1.0817^/3] 



(151) 

a result similar to Eq.(137) but without the last exponential factor. Thus the many- 
to-many amplitude would reach the unitarity limit, exp — ^ • 0, in contrast to the 
2 — > all amplitude. 

4.3. High temperature or high density transition as compared to pro- 
cesses at high energies 

As is well knownQ, high-temperature A(B+L) ^ transitions occurs classically 
(and without barrier penetration factor) in the standard model, if the temperature 
is above the energy barrier, that is if 

T > E sp (152) 

Analogously, the ground state with finite fermion density will decay without any 
suppression if the Fermi sea level is sufficiently high, i.e., of the order of the sphaleron 
massEj. 

That a similar disappearance of the tunnelling factor might occur a so in high- 
energy scattering, was the original motivation for the earliest works Ej'Ej. However, 
there is a clear physics difference here. A characteristic feature of the scattering 
processes is that the high energy of the initial channel is carried by just two energetic 
particles. In the high temperature (or high density) transitions, on the contrary, 
it is expected that the relevant initial states, at a given (high) energy and with a 
given set of quantum numbers, are mainly multiparticle states. For there are many 
more such states as compared to two particle states, hence they will dominate the 
process for purely statistical reasons. Thus although many-to-many amplitudes 



might well become unsuppressed at the sphaleron energy (see Eq.(151)), the two- 



to-mairv amplitude will remain small. A suggestion put forward by Diakonov and 
PetrovtJ'llj, that the decay rate of the state with a finite density (with the Fermi 
level /i) is proportional to A(B + L) ^ cross section (with E ~ a), seems to be 
invalidated, once this difference in the initial states are taken into accounttZI. 

5. Fermions in the valley 

In this section the whole problem will be analyzed from a somewhat different 
point of view. The behavior of the chiral fermions in the valley background is 
studied, and the A(B + L) ^ cross section is analysed by using the optical 
theorem. 
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5.1. Unitarity puzzle and Spectrum of the Dirac operator 

As already noted in Sections (1.3) and (4.1), the total A(B + L) ^ cross section 
can be computed via unitarity, i.e., as an appropriate part of the imaginary part of 
the forward elastic amplitude, 



1 



1 + 2. 



(153) 



However, a problem con cerning unitarity and chiral anomaly arises, which can be 
formulated as followsliia The optical theorem states that the cross section, 



1 



(154) 



summed over all possible X, is equal, apart from a kinematical factor, to the imag- 
inary part of the forward elastic amplitude, Eq.( |l53| ). Now consider a particular 



class of processes (154) induced by an SU(2) instanton, with 



A(£ + L)^0; A(B-L) = 0. 



(155) 



Sum over the final states satisfying (155) should give a part of the full imaginary 
part of the elastic amplitude: 



Anom Im A-r 



■X 



(156) 



see Eq.(145-146). Now, for an (anti-) instanton background, which is relevant for 
the calculation of the right hand side of Eq.(156), each right (left) handed fermion 
field has a zero mode. The standard functional integration over fermions yields a 
product of these zero modes; by going to momentum space and by applying the LSZ 
amputation one finds the S-matrix elements consistent with the instanton selection 



rule Eq.(155) 



How to calculate the left hand side of Eq.(156)? Being a part of the elastic 
amplitude, it must arise from a four point function computed in a background, 
topologically (globally) equivalent to the trivial, perturbative vacuum. To be equal 
to the total anomalous (A(_B + L) ^ 0) cross section, however, such a background 
must have a nontriviai topological structure, for instance, similar to a widely sena= 
rated instanton anti-instanton pair. More precisely, the work of Arnold and MattisE3 
suggests that one should consider something like the valley. 

The problem is that no fermion zero modes exist in the valley background!!-!, 
(see Appendix B for a sketch of the proof) T he spectrum of the Dirac operator 
in the valley background Eq.(B4|) is proven 



li 



to be the same as that of the free 
(— oo,oo) and with no normalizable 



Dirac operator, with the continuum spectrum 
zero or nonzero modes. Thus one wonders how the left hand side of Eq.(156) can 
be computed, which should somehow be dominated by the standard lefthanded or 
righthanded zero modes, in order to match the right hand side. 

Note that this " unitarity puzzle" by no means depends on the use of the fermions 
as external particles in an essential manner, although the puzzle looks much neater 
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with fcrmions. In the literature, often theories without fermions are considered, 
on the basis that the main problem lies in topologically nontrivial aspects of gauge 
field dynamics (which is quite true). The use of gauge bosons as external particles 
however does not eliminate the unitarity puzzle. For instance, the validity of the 
leading order " semiclassical" approximation (i.e. substituting the instanton solution 



in the external lines) is not at all obvious for the left hand side of Eq.(15E) (especially 
in the regime of strongly overlapping instantons), if it is more reasonable for the 
production amplitude appearing on the right hand side. 

To solve the " unitarity puzzle" , one must study the elastic amplitude starting 
from the four point function computed in an appropriate background of instanton 
anti- instanton type (which will be approximated here by the valley), and must 
show that the imaginary part of the forward elastic amplitude indeed contains 



an anomalous piece which reproduces the right hand side of Eq.(156). The work 
reviewed in the next section is a first step towards such a solution. 

5.2. Fermion Green function in the valley: widely separated instanton 
anti-instanton pair 

The four point function, 

< T^i(;r)V>2(u)$i(v)^2(v) > {A ^ ] 
V^hp rh (aOV* {y)i>2 (v) e- s /Z<- A =^ ; (157) 

S = J2 d^xi^jDiPj (158) 

3 = 1 J 

in the fixed background of Eq. (Q) , has been studied by two of the present authorsEl. 
As the functional integral factorizes in flavour the quantity of interest are (suppress- 
ing the flavour index), 

I(x,y) = J Vipvj>iP(x)i>(y)e-f dixi ^, (159) 

and 

Z = J Th/iD^eS^" 1 ^ = det D, (160) 

where it is assumed that det D is suitably regularized. 

The key point in this analysis is to introduce complete sets of orthonormal modes 

{Vn} an d {Cn }i n — 0,1,2, , for the left-handed and right-handed fermions, 

respectively. They are eigenstates of D^D^ and D^D^: 

» (e) i B) = kmCW (m = 0, 1, ...), D^Ck } = kmV^ (m = 1, 2, ...), 

= Imffl (m = 1, 2, ...), DMffl = (m = 0, 1, ...), (161) 

where 

k Q = l = 0. (162) 
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The covariant derivatives D^ a \ are defined with respect to the antiinstanton 
and instanton parts of Eq. ( |34| ) . Accordingly the zero modes are those in the regular 
gauge (for the lcfthanded mode) and in the singular gauge (for the righthanded 
one), respectively. 

The functional integration can then be defined as: 



ip{x) 



VtpVip = Yi da m db n ; 

m,n— 

oo 



oo 
m=0 



(163) 



n=0 



As is clear from the way these equations are written, the system is first put 
in a large but finite box of linear size L ( such that L 3> R, p) so that all modes 
are discrete. After the derivation of Eq.(167) below (i.e., after the sum over the 
complete sets is done), however, L can be sent to infinity without any difficulty. 

The two point function I(x, y) can be written as 



I{x,y) = 4etD(x\D- 1 \y) 

= det D {(x\a, 0)(a, 0|£» _1 |i, 0)(i, 0\y 

+ ^(a;|a,0)<a,0|5- 1 |i,n)<i,n|y) + 



(164) 



J2(x\a,m)(a,m\D- 1 \i,0)(i,0\y) 



E 



x\a, m){a, m\D 1 \i,n)(i,n\y)} 



the term proportional to the product of the zero modes has been singled out. We 
wish to compute I(x,y) at small p/R. To do this, first let us write 



D = 



( d 


Vl 


v„ 


"\ 


Wi 


111 . 


X\ n 




X m \ 


Xmn 




V : 






■J 



(165) 



The idea is that the matrix elements involving either of the zero modes, d, v n , w m , 
are all small by some overlap intergrals while the matrix elements X mn are large 
because the wave functions of non zero modes are extended to all over the spacetime. 
The inverse matrix D~ l is given by: 



(D- 



l/(d-vX~ 1 w) 

cT 1 + d~ 2 v m (X~ 1 ) mn w n 4 

(X--W&V)- 1 = X-\l 
a 



1 



-w <g) vX~ 



= (X l ) mn + d \X l ) m iWiv k {X 1 ) fe „ + ---, 



)0n — 



-d-^iD- 1 : 



In ; 



(D- 1 )^ = -(D-^ooX-lwk, 



(166) 
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where X^ 1 is the inverse of the submatrix X in the space orthogonal to the zero 
modes. 

Inserting Eq.(166) into Eq.(164) and after some algebra one finds a remarkably 
simple (and still exact) expression for I(x,y): 

I(x, y) = dot X{(x\a, 0) - (x^C^, 0)} {{%, 0\y) - (t, ^BX-^y)} 

+ ActD{x\X- 1 \y), (167) 



where C = C^oy,; B = B^a^ in Eq.(167) are defined by: 

D (valley) = rfa) + Q £,(,') + r 



(168) 



The function C M is the (modified) instanton field while B^ is the (modified) anti- 



instanton. 



Eq.(167) displays nicely the main features of the two point function in the valley 
background. The effect due to the zero modes is separated and everything else is 
expressed by the smoother two point function, 



S' x , y = {x\X- l \ 



(169) 



S' x y can be assumed to behave at large x and y (with Xi and x a fixed) as 

U\x)S F (x,y)U(y), 

&(j< ) fj. 



x,y 



U(x) 



(170) 



where S F is the free Feynman propagator. This behavior is suggested by the fact 
that the valley field has a pure gauge form at large x, 



A (valley) „ IrfdpU ~ 0{-). 

9 % 



(171) 



To proceed further one assumes that 

dct X/ det 8 = const. (172) 
as p/R — > 0. Next the ratio det-D/ detX can be estimated as follows (see ( ]166| )): 
detD 



detX 

where use was made of 



= (CD _i )oo) ~d~ const. p A /R 6 , 



(173) 



d = A>o = (», 0|C|o, 0) = / $\zyC(z) rtfXz) ~ p 2 /R 



Combining Eq.(172) and Eq.(173) gives 

det D 
det 8 ' 



p 2 /R 3 



(174) 
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With Eq.(174) and Eq.( |l70| ) in Eq.(167) one can estimate the amplitude and 
the leading contribution to its anomalous part. 

There is a subtlety here. Because the background (in the gauge of Eq.(p4f)) 



reduces asymptotically to a pure gauge form, Eq.(171), and does not vanish suffi- 
ciently fast, the standard LSZ procedure cannot be used to extract the amplitude. 
The correct procedure is to go to a more physical gauge 



^(valley) = {j^valley) + _^){/t ; 



U(x) 



yJ{x-X a ) 2 



(175) 



by appropriately transforming the fermion four point function, before the standard 
LSZ procedure is appliedti^. 



The contribution of the first term of Eq.(167) the elastic amplitude is found to 
be, to leading order in p/R, proportional to 



p 2 exp(zp • x a ) exp(-zg • x t ). 



(176) 



One finds also a correction proportional to p/R times this factor, coming from the 
terms containing the function B and C in Eq.( 167 ). These terms are precisely what 
one has been looking for: the part of the elastic amplitude, whose imaginary part 
is going to match the right hand side of the unitarity equation, Eq. (|l56| ). 



The second term of Eq.(167), in contrast, mainly gives rise to non- anomalous 
discontinuity, associated with A(B + L) = intermediate states. 

Thus at least for the valley corresponding to widely separated instanton anti- 
instanton configuration, the elastic amplitude contains the piece which reproduces 
the square of the production amplitude computed in the single instanton back- 
ground. 

The general recipe for calculating the left hand side of the " anomalous" unitarity 



relation, Eq.(156), is however not yet known to the best of our knowledge. 

We have thus far considered the behavior of fermions in the valley which corre- 
sponds to a particular kind of instanton anti-instanton pair configuration, oriented 
in the maximally attractive directions and with a particular interaction termO. 
Such a field is supposed to be of particular relevance in the problem of baryon 
number violation, as reviewed in Section 1.2. 

A related problem of the fermion propagation in a background of simple sum 
of widely separated instantons and anti-instantons (in singular gauge), was studied 
earlier. Such backgrounds mjgjlt be important in model ling some feature of the 



122 



123 



It is Lee and Bardee 



n@ 



who, in such a con- 



physical vacuum of QCDe 
text, showed for the first time that the fermion propagator can be in some sense 
dominated by a term proportional to a product of the standard zero modes (for 
widely separated instanton-anti-instantons) , in spite of the absence of any fermion 
zero modes in the background considered. 
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Also, the method used by Lee and Bardeen (they work with the equation satisfied 
by the fermion Green function, rather than doing the functional integration ) is very 
different with our own; of course, our formula reproduces theirs to the leading order 
in p/R, if one uses the simple sum of instanton and anti-instanton instead of the 
valley and sets A" 1 ~ S F + (S a - S F ) + (Si - S F ). 

5.3. Overlapping instanton anti-instanton pair 

The valley background used above reduces to the vacuum field at R — 1 . One 
wonders whether the transition to a purely perturbative field occurs gradually and 
only terminates at precisely R/p = , or it takes place abruptly at a finite value of 
R/ p, probably of order of unity. 

There is_.a strong indication that the latter possibility is realized. The first 



indicationliiil comes from the numerical and analytical study of the integrated topo- 



logical density, 

C(x 4 ) = -J 4 Jd 3 x^Tr F^F^ = Afcsfa) - Af C s(-oo), (177) 

as a function of x± for several values of R/p, for the valley background of Eq.(p4|). 
(See Fig. 11.) Ncs{%a) is the Chern Simons number 

N C s{xi) = - [ cPx-^e iijk Ti {F i} A k - ^A^.). (178) 



16tt 2 " 3" 
The instanton and anti-instanton are situated at (0, R/2) and at (0, —R/2), respec- 
tively. 

It can be seen from Fig. 11 that the topological structure is well separated and lo- 
calized at the two instanton centers only at relatively large values of R/p, R/p > 10: 
for such R/p, C[xi) reproduces locally the situation of single instanton background 
(near X4 = R/2) and that of anti-instanton (near X4 = —R/2). Vice versa, for small 
R/p < 1 the gauge field is seen to collapse to some insignificant fluctuation around 
zero, not clearly distinguishable from ordinary perturbative ones. (These state- 
ments are, admittedly, a little vague and not very precise. See the next subsection 
for an attempt to make them more precise.) 

In Fig. 12 is plotted also the behavior of the maximum of each curve, corre- 
sponding to C(0), as a function of R/p. The behavior of C(0) is powerlike both at 
large and small R: 

C(0) - 1 - 12(p/i?) 4 , R/p->\, 

C(0) ~ \(R/pf, i?/p«l. 

(Actually, the exact expression for C(x4) in terms of X4 can be found by using 
conformal transformationstij) . In particular, 

C(0H3(i-j) 2 -2(i-V (179) 

Z + 1 2+1 



'The triviality of the valley for R = is a feature independent of the choice of the weight 
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where z is defined in (|35|).) 

The behavior of C(0) supports the idea that the transition to a perturbative 
background is a sharp one. In particular, the quadratic behavior of C(0) in R/p 
found at small R/p is nothing but the reflection of the perturbative quadratic fluc- 
tuation. This is particularly transparent if one uses the gauge 



11! 



in which 



Afj, oc R. 



Altogether, Fig. 11 and Fig. 12 indicate that the instanton and antiinstanton 
start to melt at around R/p ~ 5 and go through the transition quickly, the center 
of the transition to a purely perturbativereeime being at around R/p = 1. The 
known behavior of the valley field actionE3c.3, which sharply drops from near the 
two instanton value Air/a to zero around the transition region, R/p ~ 1, (see Fig. 
3) is perfectly consistent with this conclusion. 

5.4. Level-crossing in the valley 

These discussions, though convincing intuitively, do not provide a quantitative 
answer as to the value of R/p at which such a transition occurs. The best thing to 
do would be to analyse the fermion four point function directly around R/p ~ 1, 
which however appears to be a difficult task for the time being. 

The next-best thing to do is to study the level crossing of the chiral fermion in 
the background of the valley, and see whether a qualitative change occurs at a finite 
value of R/p . The idea is the following. 

As is well known, one way to understand the anomalous chiral fermion generation 
by a topologically non trivial gauge fields, 

AN R - AN L = J a^x-^F^F^ = tf cs (-oo) - Wcs(oo), (180) 

is through the spectral flow of the eigenvalues of the Dirac Hamiltonian H, where 
7i is defined by 

i llM D„ = i l0 {D +H). (181) 

For instance, in the case of a single instanton, one finds (in the gauge Aq = 0) that 
Aj(x, — oo) and Aj(x, oo) are gauge equivalent: the spectrum of Ti is the same at 
t = ±oo. This however does not imply that individual levels Ej(t) are the same at 
t = ±oo. In fact, the existence of o ne right handed zero mode (which is known to 



be there from the index theoremlH9) , implies that precisely one right handed mod e 



crosses zero in going from t = — oo to t = oo so that (in adiabatic approximation!!!!!) , 

Mt) ~ e ~ El{oo)t ; t^ 00 ( 182 ) 

with Ei(— oo) < 0; £1(00) > 0. Otherwise, this mode would not be normalizable 
in four dimensions. 
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Furthermore the point t at which the level crossing occurs, is characterized by 
the fact that the equation Tirj = has for such t a solution, normalizable in three 
space dimensions. 

In the case of the valley background, j4j(x, — oo) and Ai(x, oo) are not only 
gauge equivalent but corresponds to the same Chern-Simons number (that is to say 
that the full integrated topological density - the Pontryagin number - is zero). Thus 
we expect that Ej(— oo) = Ej (oo) also for all j. Nevertheless, at least for widely 
separated instanton anti-instanton pair, we expect from clustering argument a non 
trivial spectral flow such as in Fig. 13a: a right handed fermion level must cross 
zero at (more or less) the instanton time position and then cross zero back near 
the anti-instanton site. In this way the physics of instanton or of anti-instanton 
(anomalous fermion generation or annihilation) would be locally reproduced by the 
valley. 

If such a picture is confirmed for large R/p, then one can ask what happens for 
small R/p. The strategy is thus to check the method of level crossing where physics 



1 2: 



explore the unknown territory, 
(see Appendix C) confirms fully 



is well understood, and then use the same too . 

The analysis made by the present authors: 
the above picture (Fig. 13a ) as long as the instanton anti-instanton distance is 
large enough: R/p > ^/a/Z ~ 1.1547.... 

(" v/4/3" is only a fit to our numerical data. However, as can be easily seen from 
Eq. ([35|) and Eq. ( |179| ) , the maximum of the variation of the Chern Simons number 
(C(0) ) reaches 1/2 precisely at R/p = \/4/3 (2 = 3). This makes one suspect that 
the critical separation is exactly \/4/3. ) 

On the other hand, for R/p less than the critical value y/4/3, no level crossing 
is found to occur. The situation is shown in Fig. 13b. If we identify the level 
crossing with the (anomalous) chiral fermion generation or annihilation (as is the 
case in the single instanton or anti-instanton background), then we arrive at the 
conclusion that the valley field with strongly overlapping instanton anti-instanton 
configurations (R/p < \/4/3 ) has nothing to do with anomaly: it is a purely 
perturbative background. The implication of this on the question of A(B + L) ^ 
cross section in the standard electroweak theory will be discussed in the next 
subsection. 

Furthermore, it was found numerically that at a crossing point t* (the value of 
t such that a normalizable solution exists), the relation 

C(t) = AMcs = \ (183) 

always holds. This result is, after all, very natural since AMcs — \ corresponds 
precisely to gauge fields sitting on top of the hill between the two adjacent vacua. 

To conclude, the results described in subsection 5.3 and 5.4 imply that the 
" instanton anti-instanton" valley configuration ceases to be topologically significant 
when the parameter R/p is equal or smaller than the critical value, y/4/3, or put 
more intuitively, when the instanton pair "overlaps" substantially. But this means 
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that the lower portion of the valley has not a well-defined physical meaning: it is 
just a sort of perturbative field, to be considered together with generic fluctuations 
around = 0. 

Such a result is not really surprising. In non-Abelian gauge theories tjifL per- 
turbative series is divergent and believed not to be even Borel summablelHS. In 
other words, the perturbative series alone does not define the theory. This is no 
reason for despair however: we know that in these theories there are also physical, 
non-perturbative effects to be taken into account. It is to be expected that only the 
sum of perturbative and non-perturbative contributions have a well-defined physical 
meaning, not each of them separately. The systematic treatment of such a sum- 
mation and a better definition of the theory, let alone the solution to the unsolved 
problem of renormalons, are however not known at present. 



5.5. High-energy A(B + L) ^ electroweak scattering 

Let us come back to the physics of A(B + L) ^ cross sections at high energies. 
Khoze and RingwaldE§ have made a simple model calculation via optical theorem, 
see Section 1.4, in which the valley solution Eq.(|34|) of the pure SU(2) Yang-Mills 
theory is used in conjunction with the Higgs contribution in the action, —irp 2 v 2 . 
By putting simply a product of the standard zero modes for the external particles 
(they actually considered a forward elastic amplitude with external gauge bosons, 
hence these zero modes are replaced by instanton or anti-instanton solution). By 
performing the integrations over the collective coordinates with the saddle point 
method, they found that the cross section grew with energy and at an energy of the 
order of the sphaleron mass, x = xkr = Sy/3/5 the exponential t'Hooft suppression 
was overcome completely. See Fig. 4a. 

Although it is a toy model calculation (the valley equation used there does not 
take into account the Higgs coupling to the gauge fields) , it is worthwhile to reflect 
upon its meaning. Since it uses the optical theorem it in principle takes into account 
the unitarity constraints: how did it manage to evade the "half-suppression" result 
mentioned earlier? 

The answer appears to be that a naive treatment of the external particles has 
led one astray. Indeed, the saddle point values of R and p (which depend on the 
energy) are such that at xkr, p = 0; R/p = 0. See Fig. 4b. The "valley" field with 
these parameters is simply a perturbative vacuum, = 0. No wonder the cross 
section is unsuppressed! 

Moreover, the analysis reviewed in the two preceding subsections strongly sug- 
gests that the transition to the purely perturbative regime occurs actually much 
earlier, when the parameter R/p reaches x/4/3. Translated into the value of the 
valley action, it means that only the "higher" portion of the valley with 

S > —£-(0.5960...) (184) 



has anything to do with the A(B + L) ^ cross section. Taking into account all 
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the other terms, we conclude that in this toy-model calculation, the A(B + L) ^ 
cross section is always suppressed at least by 

VMB+LW < e-^ ' 3185 -). (185) 



Generalizing the lesson learned above, one can give an argument for the "half 
suppression" result which does neither depend on the particular valley trajectory 
Eq.(|^) nor on its unjustified use, but only on a number of general assumptions. 
The assumptions needed are the following. 

(1) There is a semi-classical Euclidean, real background (or more properly, 
an ensemble of them as in the valley trajectory) which dominates the four point 
function, to be continued and LSZ-amputated to yield the "anomalous" imaginary 
part of the clastic amplitude. 

Although this is just an assumption - quite a nontrivial one - , and it might 
sound rather naive, after all those discussions about the complex saddle point fields 
reviewed in Section 2, there is a clear advantage here. As one uses the optical 
theorem here, all final particles are implicitly already summed over (unitarity!); the 
only effect of the external particles are four fermion (or bosonic, if one prefers to 
study purely bosonic amplitude) fields, appearing as pre-exponential factors in the 
functional integral. It is much less likely here (as compared to the calculation of 
the production amplitude, to be squared and summed over the final states) that the 
external particles affect the relevant gauge background. We recall also that in the 
approach by Diakonov and PetrovE-3 (see Section 3), the relevant fields are real and 
Euclidean (although singular). 

(2) This gauge background, to be able to give rise to the "anomalous" imaginary 
part, must have sufficiently large variations of the Chern-Simons number: more 
precisely we assume that C(t) crosses the value 1/2 at least twice (once upwards 
and once downwards, as required by Q = 0). Namely, we suppose that 

Mcs(ti) -A/bs(-oo) = -, 

Mcsfo) - Mcs(po) = \, (186) 

(ti < t2,) as a (necessary) condition for the fermion level crossings (to and back, as 
in the valley with R/ p > \/4/3 ) to occur. The action of such a gauge background 
is easily shown to be at least as large as the single instanton action, 2Tr/a. Indeed, 



since Tr(F M „ ± F^f > 0, for each x, it follows from Eq.(|186|) and Eq.(177) that 



S = - I d^xTiF: 



2 J 

> \ J_Jt J d 3 xTvF^ + ~ jTdt J d 3 xTvF^ 

ftl p 1 pOO 



1 f f 1 f°° f 27T 

" ~2 J * J d3xTrF ^F^ + g J t dt J d 3 xTvF^F^ = — . (187) 
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There are also the exponent depending on the initial energy (coming from the ex- 
ternal particles through the LSZ procedure), the part of the action due to the Higgs 
particle, as well as an integration over the collective coordinates. The saddle point 
equation however normally has a solution in which each term of the exponent is of 
the same order of magnitude. This leads to the anomalous cross section, suppressed 



by something like the square root of the t'Hooft factor, in view of Eq.(187). 

This argument is admittedly not very rigorous, and should be regarded at best 
as a tentative one. Nevertheless, it concisely takes into account the two main 
ingredients of the whole problem, unitarity and anomaly in an essential manner. 
Further efforts are welcome to see whether it can be made more rigorous. 

5.6. Note on anomaly for massive fermions 

All calculations described up to now are done for massless fermions, for sim- 
plicity. In the actual world all known fermions (except perhaps for some or all of 
the neutrinos) are massive. On the one hand, one believes that the physics at high 
energies, E 3> rrif, (rrif standing for a generic fermion mass) should be the same 
as that with m/ =0. On the other hand it is not at all obvious that anomalous 
processes such as Eq.(|j) do occur with massive fermions. Consider for instance 
the description of chiral anomaly in terms of the level crossing (recalled at the be- 
ginning of subsection 5.4). The smallest (in the absolute value) eigenvalue of the 
Hamiltonian which crosses zero once - in the case of an instanton background - has 
the magnitude which is really very close to zero: i?i(±oo) = O(j^), if L denotes the 
linear size of the spacetime volume in which our system is put. If the fermion is 
massive, the positive and negative energy levels are separated by a finite gap 2m/; 
how an infinitesimal "level crossing" (from Ei(-oo) to Ei(+oo) ) can cause a jump 
from e.g, a negative sea level to a positive one? 



The resolution of thi s pa radox (due to Krasnikov et alH3 and Anselm et al.E23; 
see also Axenides et al.liM) hinges upon the very way fermions get mass in the 
Weinberg-Salam theory (the Higgs mechanism). A fermion mass originates from 
the Yukawa interaction, 

L { y± wa = g^ L HR + h.c. (188) 



if the Higgs field has a constant part, cf> = 



Now the (constrained) instanton responsible for anomalous process Eq.(0) be- 
haves as Eq. (|l^) near the instanton center. From Eq. (|l2|) it follows that the fermion 
mass term vanishes at the instanton center, i.e, precisely where the eigenvalue of 
the Dirac Hamiltonian changes sign. For this reason the massive fermion experi- 
ences the level crossing just as a massless one. "Massive zero modes" have been 
constructed., and an analogous paradox is solved in the case of the global SU(2) 
anomaly 



6. Conclusion 
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A unified physical picture seems to emerge from various different types of anal- 
yses reviewed here. A single instanton in the Euclidean formulation describes the 
vacuum-to- vacuum tunnelling and, as such, is not directly related to the A(B+L) 
electroweak transition in TeV-energy scattering processes. The overlap of a high- 
energy two-particle state with most likely multiparticle states of the same energy 
and quantum numbers within the same well, enters as a multiplicative factor in 
the amplitude. The rest of the amplitude describes the tunnelling between states 
of adjacent "wells", whose rate rapidly grows with energy, and loses the exponen- 
tial suppression altogether as the energy exceeds the barrier height, the sphaleron 
energy. 

The first factor, which describes the mismatch between the "most favorable" 
state (for the purpose of barrier penetration) and the initial two-particle state, gets 
instead strongly suppressed as the energy grows, as has been elegantly illustrated in 
several quantum mechanical analogue problems, reported in Section 3. It is indeed 
an analogue of the Landau's semi-classical matrix elemental. 

It has at high energies the same semi-classical form, exp — — as the familiar 
tunnelling factor, although it describes a transition within the same well. Ulti- 
mately, this factor seems to be responsible for the exponential suppression of the 
A(B + L) 7^ cross section at high energies found in different analyses. The 
same factor is found also in the unitarized multi-instanton approximation (compare 
Eq.(p7|) with Eq.(H|)). 

In the instanton description, the energetic initial particles must first convert 
themselves into 0{— ) gauge-bosons, Higgs, plus the original fermions, so that these 
particles, having small energies, can be absorbed by the instanton without any 
form- factor suppression. The price paid for the conversion, ~ (a) 1 /", gives rise to 
the mismatch factoiH. 

The same Landau factor seems to be at the origin of the qualitative difference 
between the A(B + L) ^ transition at high temperature or at high fermion density 
where the transition rate becomes eventually unsuppressed, and the high energy 
scattering in which baryon number violation is most likely to remain unobservable. 

All in all, consistency of the result (at least "half-suppression" of A(B + L) ^ 
cross sections) found from a variety of different kinds of analysis (Section 3, Section 
4 and Section 5) seems to be quite remarkable. This, together with its physical 
understanding (in terms of the Landau factor as well as in terms of the general 
unitarity bound), provides us with certain degree of confidence in the result. 

To be scrupulous, the final outcome of new, semi-classical approaches (reviewed 
in Sections 2 and 3) is not known as yet. However, if these analyses eventually indi- 
cated an unsuppressed A(B + L) ^ cross sections at the sphaleron energy, it would 
come as a surprise. It would be a new (and very interesting) phenomenon, which 
has however very little to do with the original sphaleron or instanton calculations, 
and for which no hints exist for the moment. But, of course, who knows? 

Barring such a possibility, then, have we made just "MUCH ADO FOR NOTH- 
ING", after all? 
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There W^asons to believe that this is too pessimistic a viewpoint. The original 
suggestionnU'fiU that the instanton cross section is corrected by an exponentially 
growing factor at low energies is now well established. If the enhancement is so as 
to convert the 't Hooft factor into something close to its square root (an optimist 
would call it " half-enhancement" ) , 

e -ir — > e -^, (189) 

this will be insignificant in the electroweak theory, but might be rather important 
in the low-energy instanton physics in QCD. In spite of the generally accepted belief 
that the instantons play an important role in the physics of low-lying pseudoscalar 
mesons (for instance, in the solution of the "U(l) problem"E23) , quantitative under- 
standing of the role of instantons in these systems has not yet been achieved. Any 
improved understanding of physics related to instantons should be helpful. The 
same can be said, with an even stronger emphasis, with regard to the instanton 
liquid model of the QCD vacuum 



12.; 



Finally, all these issues are deeply related to the large-order b reak down of pertur- 
bation theory and connected problemsE23, especially in QCDlSEM In our opinion 



no truly significant step forward has been made after the 't Hooft 's workt^j in this 



regard. It is hoped that somewhat improved understanding of instanton physics 
achieved through the study of baryon number violation in high-energy electroweak 
interactions as reviewed here, will help clarifying these issues as well in a near future. 



Appendix A Saddle-point evaluation of multi instanton sizes 



The integral over ti, in Eq.(133T35) can be easily evaluated since it is factorized: 
it leads to the substitution, 



iEU + (^)pf P l +l /4 
9 



96tH n 2 n 2 . 

F£bF) (-1,2, ...2»-i) (A.1) 



in the exponent, Eq.(134). The final integrations over the instanton sizes have the 
form, 

(A.2) 



where 



2n 

J J ■ ■■ J \\ dpi (powers of p's) exp(I^) 



W 



-*viy+ 5 £( 



2n-l 96tt 2 „2 „2 , /r - 



(A.3) 



For the purpose of performing the above integration by the saddle point method 
one can introduce the dimensionless variables Zj instead of p's, 



( ( __j! )£ ;W r i/ 6z 5/2 ; 



(A.4) 
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in terms of which the exponent becomes 



W=(^rj)Y, (A.5) 



2n _ 2n-l 

Y ^-j24 + ( 7475) ( A - 6 ) 

z=l z=l 

The saddle point equations for Zi, though simple, do not appear to be solvable 
exactly for generic value of n. However for n sufficiently large (in practice it means 
n > 6: see below) the solution can be found easily numerically. One finds 

Zi = Z2n-i+\, (1 = 1,2, ...n), (A. 7) 

and 

z n ~ = 2~ 1/5 = 0.87055056...; n > 6, 
zi = 0.72808; z 2 = 0.85185; z 3 = 0.86817; 
z 4 = 0.87025; z 5 = 0.87051; z 6 = 0.87055 ~ (A.8) 
Inserting the above solution in Eq.( |A.6| ) one finds 

Y ~ |n - 1.06386. (A.9) 

The fact that the first and last few Zi's differ from z„ - an end point effect - is respon- 
sible for the second term of Eq. ([A.9| )) and is crucial for the subsequent discussions. 
Substituting Eq.(|A9|) into Eq.(|A.5|) one finds Eq.(|l3^). 

Appendix B Proof of absence of the fermion zero modes in the valley 

The proof of the absence of the zero modes of the Dirac operator V = ij^ (9 M — 
igA^) in the valley background is straightforward. The gauge field is given by 
Eq.(|3) which can be written as 

= -favu - Spu) F v = T]%F u a a (B.l) 

F, = \ d, \ogL{x) L(x) ee { *~ Xa) *+f {x-x l+ y)\ (B.2) 
2 [x — XiY + p z 

We consider the decomposition: 

Because the valley is a superposition of an instanton (with a righthanded zero mode) 
and an anti-istanton field, (with a lefthanded one), the only physical possibility (to 
be ruled out) is that both T> + and X>_ have a single, non degenerate, zero mode. 
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Using the explicit form of tJH), one has 

V + = -d t + io ■ V + if • F - a A f • F - F a ■ f (B.4) 
£>_ = +d t + id ■ V - if ■ F - a A f • F - F a ■ t (B.5) 

where = (t,r), 

F = rf(r,t); F = F (r,t) (B.6) 

r = Vr^, a are spin operators, and f are isospin one. 

Since T> commutes with the three dimensional total angular momentum operator 
J = — irAV+tr+f, the requirement of having a nondegenerate zero mode is satisfied 
if , ? =0. 

Let us concentrate on the left-handed mode, associated_2ath T> + below. The 
most general function with j = can be written accordinglyti3 as: 

Vja = (a 2 )j a S(r,t) - i(aa 2 )j a ■ FT(r,t) (B.7) 

The equation T> + r/ = is equivalent to: 

dtS-hd r f =0; -d r S - hd t f =0, (B.8) 

where h = h(r,t) = \ and two new functions T and S, 

T=- E L?f; S = L~iS, (B.9) 

have been introduced. The consistency condition following from the above is equiv- 
alent to V(/iVT) = (here V = (d r , dt))- Furthermore the normalization condition 
reads 

/■+oo r+OO t 

d 4 xr 2 \T\ 2 = dt dr4n—f 2 <oo: 



OO 

+ 00 p-{-OQ 



d 4 x\S\ 2 = dt dv A~kt L S < oo. (B.10) 



Note that with our choices of Xi = (—-^,0,0,0) and x a = (-^,0,0,0), L(r,t) ^ 
almost everywhere on the line r — 0, so that a necessary condition for normalization 
is f(0,t) = 0. 

One is thus led to the boundary value problem: 

V(hVT) = 0, on fi; 

f = Q, on<9J\ (B.ll) 

f e H 2 (Q), feC°(Q), 

where Q. = {(r,t) S R 2 /r > 0}. 

The operator V 2 + V/iV is an elliptic differential operator with coefficients 
analytic in £1. At this point one can use standard theorems to prove that T cannot 
have global extrema on Q. Using the continuity on f2, the boundary conditions 
Eq.(B.lC) and the condition T — > 0, r, t — > oo, it is easy to prove that T = 0. 



R. Guida, K. Konishi and N. Magnoli 53 



It follows that S must be a constant, i.e, 

3 

S = const. L(x)~z . 

The normalization condition however forces S = 0. Q.E.D. 

Having disposed of the zero modes, one might then ask if the zero mode of T> + in 
the antistanton field, and that of 2?_ in the istanton background, combine somehow 
in the valley field so as to form a pair of non zero eigenvalue ±A of T>. A theorem 



by Ikebe and UchiyamaliiS however excludes this possibility. 



Appendix C Level-crossing in the valley 

The Dirac Hamiltonian TL is defined by: ij^D^ — 170 (-Do 
Decomposing 

' H+ 
H- 



H = 



and passing to a noncovariant formalism, one finds: 



H). 



H+ = 



-H_ = H, 



(C.l) 



with 



H(t) = +ia ■ V - a A f ■ F - F a ■ f 



(C.2) 



where the notation is the same as in Appendix B, Eq.(B.G). 

To simplify the problem, one can just look for zero modes of Ti.(t; R) without 
studying the detailed behavior of levels in terms of parameters t, R. 

The condition ( p.l[ ) implies that zero modes always appear in pairs of opposite 
chirality. We assume that for each chirality zero modes are non degenerate, so 
that we must look for singlets of the total spatial angular momentum (which is a 



symmetry of Tl). The most general form of the singlet 77 is given in ( B.7 ) 
The equation Hrj = then reads: 

SF S -3T- rVT + 2FrT = 0, 
-F r 2 T + r\7S + 2FS = 0. 



(C.3) 



Making the substitution 



we get a simplified system: 



S=jS, T=^f, 
L r A 



f'(r;t)=3F ^S; S'(r;t) = F ^f, 



(C.4) 



(C.5) 



where the prime means differentiation with respect to r and the time t appears as 
a parameter. 
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Normalizability of the solution implies: 

foo poo 2 foo poo r2 

oo> / drr 2 \S\ 2 = / ^-^l^l 2 ; oo > / drr A \T\ 2 = / dr— \f\ 2 . (C.6) 
Jo J o L J J r 

We wish to know for which range of the parameter R and for which values of 



t - if any - (the other parameter p just fixes the scale) the system Eq.(C.5) has a 
normalizable solution. Normalizability enforces the following initial condition (at a 
given t ) 

5(0) = 1; f(Q) = 0; if t ? y + | (C.7) 
5(0) =0; f(0) = l; if t = y+|, (C.8) 

where y is defined (|35|). 

The problem turns out to be too hard to be treated analytically: it suffices to 
note that the system is equivalent to a second orde r linear equation for T with 
12 regular fuchsian points! In passing we recal£H3 that for the case of a pure 
anti-instanton (or an instanton) a normalizable solution can be found analytically. 



Lacking a better method a numerical method was used to solve Eq.(C.5). The 
strategy adopted is as follows. 

(i) First solve the system by a power series in r at r ~ and r — > oo. It turns out 
that in each region one and only one of the two independent solutions is compatible 
with the normalizability. 

(ii) Choose, for a given t and R/p, the normalizable solution for T,S near t = 0. 



Make them evolve according to Eq.(C.5) up to a very large but fixed value of r, r\. 

(iii) Plot T(ta) as a function of t, for a fixed R/p. See whether T(r&) crosses zero 
as t is varied. If it does at some t, it means by continuity that for t very close to 
that value there is a normalizable solution. The power series solution of the coupled 
equation assures that whenever T is normalizable, so is S. 

(iv) Repeat the same procedure for different values of Rj p. 

Note that this method is somewhat similar to that used in th e pr oof of the 
so-called oscillation theorem in one-dimensional quantum mechanicsEHj. It is also 
reminiscent of the "shooting method" used by mathematicians for solving certain 
differential equations. 

The result cited earlier (Fig. 13a and Fig. 13b ) was obtained this way: for all 
values of R/p above a critical value, 

^ ~ 1.15470.... ~ y|, (C.9) 

there are two values of t for which a normalizable solution of Hi] = exists. Below 
the critical separation, no level crossing occurs. 
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Figure captions 

Fig 1. Graphs representing: a) the leading approximation, b) a soft-soft correction, 
c) a hard-soft corrections, d) a hard-hard correction, in a theory with fermions 
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(in an instanton background). 

Fig 2. The leading contribution to the exponent in the i?-term approach: double 
residue of W classical field integrated over the phase space. 

Fig 3. The action of the instanton anti-instanton valley (in pure gauge theory), 
plotted in terms of R/p. 

Fig 4. a) The saddle point exponent, b) the saddle point ratio R/p, versus the 
parameter x = E/E sp in the Khozc-Ringwald model. 

Fig 5. The lowest hard-hard correction in a pure gauge theory. 

Fig 6. Example of "squared tree" multiloop graph in the instanton background 
(blobs means tree graphs). 

Fig 7. The relevant complex-time path in the multiparticle approach. 

Fig 8. The behavior of the singular solution used in the W.K.B. approach (B and 
C are to be intended located at — oo). 

Fig 9. The contours 71 and 72 in the complex q plane, used in the W.K.B. approach. 

Fig 10. A multiinstanton chain. The solid lines indicate exchange of many bosons, 
the lines with an arrow representing fermions. For simplicity of drawing Np 
is taken to be 2 here. 

Fig 11. C(x 4 ) versus x^j p for ^ = 10 (outmost curve), = 5, = 2 (middle), 
~ = 1 an d ^ = 0.5 (innermost curve). 

Fig 12. The maximum of C(x±), C(0), as a function of ^. 

Fig 13. Schematic behavior of the level crossing for the Dirac Hamiltonian in the 
valley background, a) for R > R c , b) for R < R c . 
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